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CONSERVATION LAWS AND FERN-EQUIVALENCE IN GENERAL 
КЕШ 


By Jerzy Rayskı 
Institute of Physics of the Jagellonian University, Kraków. 
(Received January 24, 1961) 


А generalized parallelism (called fern-equivalence), fern-equivalent tetrads, and quasi- 
cartesian coordinates are defined in riemannian geometry with the help of one-parametric families 
of extremal (minimal) hypersurfaces. The same ideas are applied to secure conservation laws and 
localization of energy, momentum, and angular momentum in general relativity. 


1. Introduction 


In a paper entitled “Riemann-Geometrie mit Aufrechterhaltung des Begriffes des 
Fernparallelismus" by Einstein (1928) there was expressed the following opinion: „Die 
Riemannsche Geometrie ist dadurch charakterisiert, dass... die Ветаве zweier Linienele- 
mente, welche den infinitesimalen Umgebungen zweier endlich voneinander entfernter 
Punkte Р und О angehören, miteinander vergleichbar sind. Dagegen fehlt der Begriff 
der Parallelität solcher zweier Linienelemente; der Richtungsbegriff existiert nicht für 
das Endliche“. 

This opinion, or rather its second part, shared generally by physieists, is based upon 
the fact that, according to the well known definition of parallel displacement, the result 
of the parallel displacement is not unique but depends upon the path of the displacement. 
In particular, a vector displaced parallelly from the point P along a closed path back to 
the same point is no more the same vector. However, the above quoted opinion is misleading. 
Einstein himself denies it in the same article by looking for a definition of fern -parallelism. 

Led by his physical intuition, Einstein regarded the lack of the concept of fern-parallelism 
as a serious drawback hampering the physical interpretation and tried to supplement the 
riemannian geometry by a prescription permitting to decide which one of the variety of 
directions at a point Q is to be regarded аз “parallel” to a given direction at the point P. 

To achieve this end he introduced a tetrad field hÀ (greek indices or numerals 1,...4 
denote components while latin indices from the beginning of the alphabet or numerals 
in parenthesis (1),...(4) denote vectors constituting the tetrad). Then he defined in terms 
of h^ some analogue of the affine connection permitting to formulate a unitary theory of 


(509) 


510 


gravitation and electromagnetism. Which of the directions at Q is to be called fern-parallel | 
to a given direction at P, depends — according to Einstein — upon the actual state of the: 
gravitational and the electromagnetic fields. 

One of the aims of the present paper is to show that, besides the einsteinian concept 
of fern-parallelism, there exists another, by far simpler alternative. It is possible to generalize : 
the notion of parallelism in such a straightforward and natural way that (not only line ele- - 
ments but also) directions at any two points become comparable within the framework of | 
riemannian geometry. 

In order to distinguish our concept of fern-parallelism from that of Einstein we shall. 
speak about "fern-equivalence" and about “‘fern-equivalent directions". 


2. Extremal hypersurfaces and fern-equivalent directions 


Let a unit vector e^(P) represent a given direction at a point P in an n-dimensional 
riemannian space. It can be either a time-like or a space-like (but not a light-like) vector. 
In order to define unambiguously which of the directions at another point Q is equivalent 
(“parallel”) to that given by e*(P) let us make the following construction. We introduce 
an (n—1)-dimensional minimal (extremal) hypersurface with the following properties: (1) the 
point P is lying on this hypersurface, (ii) this hypersurface is orthogonal to e*(P). If the point 
Q is lying on the same surface then the direction orthogonal to it at the point Q is to be re- 
garded as equivalent, i.e. parallel in the generalized sense of this word, to that represented 
by e^(P). Introducing the unit vector e^(Q) orthogonal to the above mentioned hypersurface 
and directed so that &^(Q) goes continuously into e^(P) for Ρ-»0 (i.e. pointing to the same 
side of the hypersurface), the vectors e4(P) and e^(Q) can be called fern-equivalent. 

In the limit of flat space the minimal hypersurface goes over into a hyperplane and the 
directions e^(Q) and e (P) become parallel so that, indeed, it means a natural generalization 
of the notion of parallelism. 

If the point Q is not lying on the above mentioned hypersurface then we make the fol- 
lowing construction. We introduce a one-parametric family of extremal hypersurfaces 


φία”) =p, (1) 
where @ is a scalar function and p is a parameter, so that 
P(x") = Po (0) 


denotes the particular hypersurface introduced formerly. The function g is restricted ex 
definitione by the condition 


0<|h*h,|<M (2) 
where 
h, =9,ф (3) 


for every point in the whole space. The family of extremal hypersurfaces satisfying (2) 


okt 


and (3) will be called "uniform"!. Now, the versot of the vector A" 
e" = |h"h,|- "^h" (4) 


taken at an arbitrary point Q represents a fern-equivalent direction to that represented 
by the original e^(P). In this way we defined a field of fern-equivalent unit vectors ο, 

In the limit of flat space the one-parametric family of extremal hypersurfaces goes 
over into a one-parametric family of hyperplanes. In particular, a one-parametric family 
satisfying (2) goes over into a family of parallel hyperplanes. Hence, the vector field е^ goes 
over into a constant vector field. Thus, our definition of fern-equivalence is a natural gene- 
ralization of the concept of parallelism for the case of riemannian spaces. This generalization 
is as natural as e.g. the generalization of straight lines to geodesic lines. Both generalizations 
are of a similar type since both take advantage of the extremal properties inherent in the 
riemannian geometry. Geodesic lines are one-dimensional extremal manifolds whereas we 
have exploited now the (n — 1)-dimensional extremal manifolds. 

Аза by-product we get another generalization of the concept of straight lines, quite 
different from geodesic lines: the orthogonal trajectories to a one-parametric family of 
(n—1)-dimensional, extremal, uniform hypersurfaces. 

According to a theorem of Bianchi (1903) and Bompiani (1921) a necessary and sufficient 
condition that the family of hypersurfaces (1) be extremal is that their orthogonal trajec- 
tories determine a correspondence between them which preserves volume. From this theorem 
it follows that 


ο == 0. (5) 
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Thus, our family of extremal hypersurfaces is characterized by a curl-free vector field (3) 
whose versors form a source-free field. In the next sections we shall exploit mainly the 
analytical properties (3), (4), and (5). 


3. Fern-equivalent tetrads and quasi-cartesian coordinates 
Let us introduce four independent uniform families of extremal hypersurfaces 
g^ =p, απ)... (4), (6) 


in the four-dimensional riemannian manifold to be interpreted physically as space-time, 
and assume three of the families to be time-like (for a=(1), (2), and (3)), and the fourth to 
be space-like (for a=(4)). Introducing the vectors hî and their versors ей according to (2), 
(3), and (4), we get а source-free tetrad field ` 


e =0 (5) 
satisfying the conditions of normalization but not of orthogonality 
ead = y (7 


1 Uniform families of minimal hypersurfaces surely exist in spaces which become flat for points tending 


to infinity in space-like directions. 


but 
y? = y £ 0 for a x b. (73 


Such a field of fern-equivalent tetrads will be useful for physical interpretation. 

Two observers situated at P and О whose world lines are just tangential to two orthogonal 
trajectories of a uniform family of extremal hypersurfaces p® = р, i.e. whose directions are 
(momentarily) е*®(Р) and e“®(Q), may be regarded as being at rest one relatively to the 
other. Thus, it is possible to introduce the concept of a quasi-inertial frame of reference 
by forcing the bodies of reference to move along the orthogonal trajectories of a uniform 
family of space-like extremal hypersurfaces. A transition from one quasi-inertial frame of 
reference to another is done by going over to another family Фа = р of Ще same type. 

It is also possible to introduce а quasi-cartesian system of reference by defining the 
invariant coordinates 

Ρ 


ΤΡ = fe WE ds. (8) 


$ 
Р, 


Owing to (3) the above integrals are independent of the path of integration and сап be 


written also as 
P 


“ = | аре = 9°(P) — φρο. (9) 
Po 
If the space-time is of such a type that for Р-> со in space-like directions the manifold be- 
comes flat and if we assume additionaly 


ηλ. > Nab (10) 


for P> in space-like directions, then the coordinates r^ become a natural generalization ` 
of cartesian coordinates. A transition between two quasi-cartesian systems of reference 
r,>r,is achieved by m: iking a Lorentz transformation of the basic vectors at space-like in- 
fu. 


Let us introduce the reciprocal tetrad field ey, by 


A 


езер = дв whence eie as OF, (11) 
It satisfies the relations 

wl ere Gu а δα 19)‏ 5700 ۸م 

{coe γα У Уд = ος: ( 


Differentiating the second formula (11) with respect to a“ and taking into account (5) we 
get the useful formula 


е 9 ея = Г, (13) 
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4. Local conservation laws in general relativity 


In spite of numerous efforts in the past forty five years the problem of local conservation 
laws for energy-momentum and angular momentum in general relativity has not been 
solved yet. In view of this situation several physicists represent the opinion that local con- 
servation laws cannot exist since the riemannian space-time is neither homogeneous nor 
isotropic. We are not convinced by this argument. Conservation laws are violated if the 
| lagrangian depends explicitely on the arguments х“, which is not the case in general relativity. 
"The lagrangian of general relativity depends only implicitely upon the arguments a" through 
the intermediary of the field quantities g,,. Conservation laws in special relativity are deriv- 
able either directly, by guessing the form of the conserved quantities, or with the help 
of the well known procedure of infinitesimal rigid displacements. Such rigid displacements 
are believed to be impossible in the riemannian space. In fact, it is impossible to displace 


rigidly the riemannian space into itself, i.e. without going beyond it, but surely it is possible 
to displace it rigidly within the ten-dimensional embedding space. If all material bodies 
(fields) including the gravitational field. itself are rigidly displaced in the embedding space 
then the observers can use new coordinates and new tetrads displaced so that no change 
will be noticed except for an immaterial change of поѓалоп. Therefore the existence of local 
conservation laws is to be expected also within the framework of general relativity. 

In order to show the possibility of securing localization and conservation of energy, mo- 
mentum, and of angular momentum we shall not take advantage of the possibility of 
performing rigid displacements in the embedding space, but we shall use the tetrad 
technique and the concept of fernequivalence. 

From the purely technical point of view the trouble with the conservation laws in 
general relativity is that the covariant continuity equations 


Ge = 0, (14) 


1 
where — G” is the symmetrical energy-momentum-stress tensor 
a ; j 


Ge = - | Re — Larn) ss ee 


are not suitable for an application of the Gauss-Ostrogradzki theorem. However, as is well 
known, there is no difficulty whatsoever with the conservation and localization of the electric 
charge in general relativity in. spite of the fact that in this case we have also to do with 
a covariant form of the continuity equation. The technical reason for the lack of any 
difficulty with the charge and an apparently insurmountable difficulty with the energy- 
-momentum is due to the fact that covariant derivatives mean something else for each case 
of tensors of different ranks. The case of tensors of the first rank Т^ is, so to say, natural 


for deriving conservation laws from covariant continuity equations, 


о (16) 
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Indeed by integrating (16) over an arbitrary four-dimensional volume О, we get an expression | 
to which the Gauss theorem is directly applicable 


di f 6,77 ТОТ V — gdx = J 3.7” V—g)dx =>; ф T^y —gas"e (17). 
2 2 


(uveo) $ 
where use was made of the well known formula 
әү —&g-l,. (18) 
Introducing the element 45, dual to (75% 
dS, — / — gds", (19) | 


splitting the closed hypersurface into two space-like hypersurfaces Y, and 2, changing the 
direction of the normal on У, we get 


|5 = [ Tras (20) 
whence the quantity 
T= | T*dS, = if T"dS (21) 


is independent of the choice of the space-like hypersurface, i.e. is conserved. The normal 
component 7” is interpretable as the density of T at the points of X. 

In the above derivation the form of the covariant derivative of the vector was explicitely 
taken advantage of. In the case of tensors of higher ranks 7", where the covariant deriva- 
tives involve further terms with Christoffel symbols, the analogous derivation obviously 
breaks down and it ıs impossible to infer the conservation laws. 

The above discussion not only shows where technical difficulties arise from but, at 
the same time, suggest the way out of the difficulties: The tensor of energy-momentum 
is to be replaced by a suitable set of vectors. This end can be achieved by using the tetrad 
technique. By means of projections upon the legs of the four-leg (tetrad) A?" we can replace 
the energy-momentum-stress tensor E" of the complete system including substantialistic 
matter as well as gravitation by a set of four fourvectors 


Wee — Һан ри, (22) 


/ 


These vectors may be called the energy-momentum "currents", 
If we succeed in satisfying continuity equations of the following form 
WU =0 (23) 
then, repeating exactly the derivation (16)—(21), we shall find strictly conserved and localiz- 

able quantities 
‘a j an 
W*— | w as, (24) 
S 


to be interpreted as the total energy-momentum components with respect to the tetrad field 
ai 2 2 : E 
Л“. The normal components W а will be interpretable as the corresponding densities. 
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Thus, in order to guarantee conservation and strict localization of energy-momentum 
in general relativity, we have to construct the "currents" satisfying the continuity equations 
(23). It is obviously possible to find a symmetrical tensor W , satisfying (23), i.e. to guarantee 
local conservation laws for energy-momentum. However, as there are ten unknown functions 
W „and only four equations (23), the problem is still underdetermined. In order to determine 
W,, from (23) some more information about the structure of IP is needed. 


WT can be defined as 
Wy = № 0 hi (25) 


where OF is the generator of infinitesimal translations of quasi-inertial coordinates г“. 
The generator @ is identical with the well known affine tensor of Einstein if (and 
only if) the latter is expressed in quasi-inertial coordinates. But 0% has а covariant 
meaning: it is a set of scalars with respect to general coordinate transformations and 
transforms like a tensor with respect to transformations replacing one system of families 
of extremal hypersurfaces by another system of extremal hypersurfaces. 

The tensor Л“ can be symmetrized by the well known procedure of Belinfante 
which secures local conservation laws of the quantities | 


Mae ер bu ра we 


representing angular momentum. 
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THE GENERALIZATION OF DIRAC'S EQUATION II 


By Jerzy LUKIERSKI 
Institute of Theoretical Physics, University of Wroclaw 


(Received January 23, 1961) 


The quantized Lagrangian formalism, which leads to the generalized Dirac equation in- 
troduced in I, is investigated. The theory is invariant with respect to the 12-parameter € x €" 
group. The Lagrangian is described by means of the two real isovectors: the mass-isovector 
Ж. (introduced in 1) and a second isovector EES not occuring in the equation. The isovector er 
is necessary for quantization invariant under СХ ©’. It is shown that the Jauch field may be 
obtained and generalized after a special choice of Esu: 


Introduction 


In the present paper we continue some considerations published recently (Lukierski 
1960 !) was to introduce the generalized Dirac equation, invariant under to the 12-parameter 
group C x С’. We have stated on a generalized Dirac equation described with help of an 
isovector x,, standing for the scalar mass in the usual Dirac equation. It is also shown that 
if the mass-isovector is chosen as follows (see (1.1.15)): 


4. = (0, 0, — m, 0) (0.1) 


зи 


we obtain the conventional Dirac equation. 

In the following considerations we go over to the quantized theory. The most general 
Lagrangian, invariant under the C x С’ group, bilinear in the field operators ` ae: 
Yj... and linear in the derivatives (Sect. 1, 2) leads to the equation investigated in I (Sect 3). 
The formulation of the Lagrangian formalism with help of the 8-spinor notation yields 
a possibility of quantization (Sect. 4,5). 

The Lagrangian for our general equation is not determined uniquely. We obtain three 
new degrees of freedom, described by means of the coordinates of a second isovector €, З 


which satisfy the following relation (Sect. 7): 


ee 0 (0.2) 


1 Quoted hereafter as Г. 
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The introduction of the second isovector Ε΄, permits to interpret the Jauch field: 
(Jauch 1954); this field is obtained as a special case of our generalized field (Sect. 6). Schwin- 
ger's choice of the Lagrangian (Schwinger 1948) implies (besides (0.1)): 


£3 ο ος (0.3) ) 
and for the Jauch field we get 
є. = (В, 0,0, х) αἲ--βὂ--] (0.4) ) 
We see also that the Jauch field for which 
(v(x), v(x)) = — 1558 (х — x), 
р(х), p(x’)} = — 105, (х — a"), (0.5) | 
where ф = yy, and 
О<о<1 о=о* (0.5a) | 


can be generalized to the 2-parameter Jauch field, for which instead of (0.5a) 
ог ag" (0.5b) 


Let us finally enumerate some fields introduced in the last years as generalized Dirac 
fields. These fields are described with help of our generalized equation: 
A. The Ozaki field (Ozaki 1960), introduced by the following equation 


(y,9, — em} yp = 0 (0.6) 
where f = В*, can be obtained by putting 
κι, = (0, 0, cosh B, sinh В) (0.7) 


Equation (0.6) was also introduced in another form by Woythousen (1958). 
B. The Lurie field (Lurie 1959), obtained from the conventional Dirac field by using 
the following 4-parameter subgroup of С” 


‚Bys 


vw = h” (ay + Буур) (0.8) 


where 
ΑΛ" |4 pp 1 θεα 
and the upper index “c” denotes charge conjugation. 
C. The Jauch field (Jauch 1954), described by means of (0.1) and (0.4). 


D. If, according to a proposal of Rzewuski (1958), we introduce in spinor space the 
following 1-parameter group 


a: =e (0.9) 


ία" asp 
we obtain the Nishijima field in the conventional 4-spinor Dirac equation, 
(γι, -- ет} =0 (ëss oft (0.10) 
a field introduced for leptons Бу Nishijima (1957). 
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1. Bilinear Forms Constructed. from Spinor-Isospinors 


It is easily seen that from the simplest true non-hermitian representation vices CO the 
o H 
С x C' group (see I) one can construct? only the following five bilinear tensor forms in € x С”: 


== AU (1.1) 
απο s EE (1.2) 
Es = ος (1.3) 
με = e (1.4) 
ноут — И. (1.5) 
where 
и = Ate, E E ο ο ο) 
ыы οι) (1.6) 
Ze - AIR ZEN eal (1.7) 
с =o — metric tensor in spinor space (Corson 1954), σ,., σ., --- well-known 2 х 2 Pauli 
matrices. 
It may be proved that in classical theory 
fire acc ES νετ (1.8) 
and 


TT? sU PT 


T x ГО. (1.9) 
From (1.9) we obtain Takabayashi’s Lorentz frame of }-spin particle (Takabayashi 
1958a, b). 
Let us consider the non-commuting field operators Hl, У. Because 
Uer = Eugen Ee, (1.10) 


and analogously for 7',,*,, the antisymmetrical tensors T,,,, T,,, and Т, are fully 


рт oy? 


determined by their space and isospace components: 


Tied ena (1.114) 
Ta; = i i Г (1.11b) 
T ias = 2 η L (1.11c) 


? We indicate the components in € with an index before the semicolon, and the components in €' — 
with an index behind the semicolon. Because all the spin-tensors investigated below belong to usual or iso- 


-space, it is sufficient to indicate with semicolon the tensor indices only. 
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We see also that the transposition operation (classically equal to the identity transforma4 


tion) leads to the following equalities: 
Тее Ты = Tou 
Te Т. T,=-T, (1.12) 


21] 1); Us 


T Тау Т pô 
lap ies Ом Py; 


из» ER 


where only T,., and CA, are independent. 
Taking into account the hermiticity of the tensors (1.1-5), it is easily seen that only; 
T,,,is hermitian. Thus from each form (1.1-5) we obtain two different hermitian tensors,: 


which are scalar or pseudoscalar with respect to the transposition operation (1.12): 


peu SL Tu TO=Z(T-TH (1.132) 
i | E 1 Tes 

ECH MEX (Tig — Τη) Тул? к, (Т, + T) (1.130) 

1 
T = = fus η Ta) (Ee Es 9 (Tix. qu Τη) (1.130) 
TAS) = $ (Tap TL) (1.13d) 

Ti, = Ke t. ik; st Te я) 

І T+ 
Th Se 2 (Т.я = Т; ῃ)- (1.13е) 
Let us introduce the following transformation 

KY, om UE REH (1.14) 


which commute with the unitary subgroup U x U’. If we put 


Kn (1.15) 


(А . pat . | 
we see that T describe tensors, and n ) — pseudotensors with respect to the: 
inversions of time and isotime, defined with help of formula (1.14). 


2. The General Lagrangian 


The most general bilinear Lagrangian leading to the first. order equation for Ч, 
which is invariant with respect to the group € x С”, is the following one: 


C= Li+ Lat L (2.1) 


where 
of = 


> (+ BP, Er) | (2.2) 
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а 

ABC ri V pP 68 + H.C. (2.3) 
b . el 

dE >= SE a; p ο ai ep HEO: (2.4) 

where «, В are real, and a, b complex. 

If we introduce N 
6:89 = ο 6 (2.5) 
after putting 

2388 — (gi (2.6) 


we have from the hermiticity condition for £4. 


ee елт 6. (2.7) 
If we suppose that 

p pe (2.8) 

the tensor E,,,, obtained from Ей, with help of the formula 

По 801 В 
Е δ: 2 I ὃ τς (2.9) 
may be assumed to be antisymmetrical and real?: 

E m= Ey Hp Eu (2.10)' 


Further limitations on (2.1) follow from the consistency of the two equations 


Б ifi ers gif Wt, + ар «if + b ET e = 0 (2.11) 
OF asp Ё 
д AP | 
L Ξ-ἰαδὰν; RP, + а* (DENE b* LE dr (si? = (2.12) 
dÉ, 


where the variations are taken from the left. Writing out the complex conjugate equation 
οἱ (2.12): 
Ισ APR DER apes — 0 (2.13) 


we obtain following the restrictions on the parameters mos 
Neap 0 00 
В) α--β a=0 b#0 


(2.14) 


3 It is also possible to assume the self-duality condition, which implies three independent complex com- 


ponents satisfying (1.11b). 
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which lead to the following two Lagrangians 


te = 5 EMA αν ах у 69°" аря) dE E Lue A ues = нс.) (2.15) 


У . . . р . b ч А 
Ly = (евр, — Ψ, PL) + (s EP p4 μ.ο) (2.16) 


Using (1.13-15), И is easily seen that if 
a=-a* b=-b* (2.17a) 
we have 
Kf phot C. 
πρ i-r. (2.17b) 


thus the Lagrangians L4, „Св have the determined parities. 
The equations obtained from (2.15) and (2.16) are 


Вау, = ibep +, — aU ۴ذ“‎ ϱ (2.18) 
REMY, ο +, + HEP PER — 0 (2.19) 
If we put 

€^ 2.2 yet, (2.20) 

we have 
Re — E? Re (2.21) 

and if 

Ts = AEN #0 (2.22) 


the equation (2.19) can be obtained from (2.18). 
Introducing (0.9), where 


[с, a] = le, a] =0 (2.23) 
we see that under this transformation in (2.1) 
a>a = e ۾‎ 
ὃ ج‎ b! = ε--αῦ, (2.24) 


and taking into account the group @, the conditions (2.17a) might be always satisfied. There- 
fore, the equations (2.18-19) are invariant under (1,13), 


3. The Correspondence with the Klein-Gordon Equation 


Let us assume that 


From (2.18) we obtain 
En i9, e.g АШ Sa в τῳ 


and, using 


we have from (2.18) and (3.2) 
αρα ος 0 


α 


The equations (3.1) and (3.4) are equivalent, if 


and 


It is easily seen that after multiplication of (2.18) by ie,;,, we obtain 


ergëtt. — Ay p PP = 0 


where 
A= ας 
[αι 
and 
us 


accordingly to the results of I. 
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(3.1) 


(3.2) 


(9.5) 


(9.4) 


(9.5) 


(3.6) 


(8.1) 


(3.8) 


(3.9) 


(3.92) 


(3.9b) 


Let us investigate the correspondence of (2.19) with the Klein-Gordon equation. If 


Ej A (E) 


(3.10) 
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the complex conjugate equation to (2.19) is 


— i9, c éd — b E GPS = 0 (3.11) 
Now we multiply (2.19) by ЗА and (3.11) by — ibE ^ and obtain by addition 
ie’? OWA + БО; (eg Ea — EP? + ib*EPESVS-O0 (3.12) 

From the relation (2.9) we: have 
EJES =r (REFE E (3.13) 

where 

be 3.14 
s = GH Euvix E six (: . ) 


By virtue of (3.12-14) we see that if 
es E, Ed (3.15) 
and 
Rb = 0 (3.16) 
equation (3.12) is equivalent to (3.1). In this case 


2 4 
mi BPE sur Eye (3.17) 


4“ 


€ 


The supplementary condition (3.15) may be described by help of the isotensor compo- 
nents €, E,,,. Because (Corson 1954) 


SN Piet E Ka T Єл»), „сд (3.18а) 
Te Oaks = Ee Dé баіо (9.189) 
where 
ο = дд, >= да, (3.19) 
we have 
eji E^, = e, (E, iE, „до, (3.202) 
Ke ES e (£s, xd IE и») 6,18 (3.20b) 
and (3.15) implies that 
ep e (3.21) 


It is shown in Appendix I that the limitation (3.16) follows from (3.21). Therefore the 
supplementary condition (3.21) describes all limit 


in. 


alions on the isotensor components 


525 


Let us multiply (2.19) by is We obtain thus equation (3.9), where 


η (3.224) 
and 
ра 
A=-i Л (3.22b) 
Considering (3.204) and (3.222), we see that our supplementary condition (3.21) is equiva- 
lent to the reality condition of the mass-isovector κ... It may be also noted that in the equation 


(2.21) the isovector e,, from R, can be chosen real only if relation (3.21) is satisfied. Indeed, 
from (2.20) we have 
a 


д = 37 26/385 
7 BEA, ER, f° X (3.23) 


or, by virtue of (3.10) and (3.20b) 


є, CE E) (3.24) 


= a 
» — БЕЯ, BB, 


By assuming (3.21), we obtain 


ipee с = ep (3.25) 


where if 
Im ig E (3.25) 


the isovector e,, is also real. 
Finally we conclude that £, and Lp with the supplementary condition (3.21) implies 
the generalized Dirac equation introduced in I. 


4. The 8-spinor Notation 


Let us introduce the 8-spinor Φ (see 1.4.6), which satisfies the following supplementary 


condition 


D=irl; LO = АФ = АФ (4.1) 
where 
qug a ar 
i 1 4 ш 2 γι (4.3) 
yi = 0223 ασ XT 
and 
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Because 


Ф = T, = 2 (P, — реу) (4.4) 


we obtain from Sect. 1 
TCD Фо ΤΟ = iDT,D 
Te? = +6 Т? = = Фт, Ф 


TOO Aen ft) = Фе Г.Ф (4.5) 
TD = = iT. Tyr, ð 
Dez "ФГ v, o 
where 
т, = (т;, — Г) (4.6) 


The relations (4.5) may be written in an arbitrary y,-representation. 
If we substitute in (2.15-16) the tensor forms (4.5), the Lagrangians may be written 
as follows: 


£, = DES DC EE Φ exp (1*Г.) Ф, (4.7) 
Lp = ФГЬГ,Ә„т„є„Ф + |b| ® exp (ВГ) @Г„Е-+ D) 20, (4.8) 
where 
a=|ale* b= ей _ (4.9) 
EL DIE Е со О 
а Бы (4.11) 


D= (Ei. E og. Е 93) 


The equations obtained from A rA and Zo are as follows: 


Rı% = (Гл, + |a| =) 6 — (4.12) 
Ra Ф = {Г.Г лье, + |b] οὔτι (Г. SS 12-0 (4.13) 

From (4.1) we have 
[A, Ву] = (4, R,} = 0 (4.14) 


The equations (4.12) and (4.13) are identical, if (assuming œ = f) 
А Жылы е RET 
Е, = e GP E-- D)T- В, (4.15) 


The 8-spinor formalism is very useful in investigations of the commutation relations, 
as the field equations (4.12-13) are linear. 
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5. The Commutation Relations 


The conventional method of Jordan-Pauli-Fierz (Jordan and Pauli 1928, Fierz 1939, 
. Fierz and Pauli 1939) for quantization of a free field is based on the definition of canonical 
variables, We shall quantize our field without the canonical formalism, using methods, 
developed by Umezawa and Takahashi (Takahashi and Umezawa 1953, Umezawa 1956). 
In this method we obtain the commutation relations directly from the field equation, if 
this equation can be reduced to the Klein-Gordon equation. 

Let us have (see Takahashi and Umezawa 1953, eq. (2.1)) 


В,,(9,) 0, (9) =0 et 


where (5.1) is obtained from the Lagrangian by taking the variation 6Q, from the left. 


Thus if 


Rag (9,) Dz, @,) = (O — т?) д,в (5.2) 
we obtain (see Takahashi and Umezawa 1953, ©. (Z) 
[ζω Qala =19,,(8,) A (x — 2’) (5.3) 
In our case 
Q=] 0,29, (5.4) 


If we assume, in accordance with (3.6), 


і ем (5.54) : 


and 
HS suv 61 eo — 62 (5.5b) 
where n St Lt 
we get from (3.5) and (3.17) 
m = |а] (5.6a) 
and 
m = eies |bl, (5.6b) 
which implies 
6162 = 1 (5.7a) 
or 
61 = Ca = (5 (5.7Ъ) 


From (5.2) and (4.12-13) we have (assuming (5.6)) 


D, = Гал, — me (5.8) 


ev’ sv 


D, = — ell el дт SE те- ГГ. E+ D) τὶ (5.9) 


HI sp 
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It is easy to prove that 


If we define 
K gp (x — х) = (4 (x), Pg (^) 
from (5.3) and (5.10), we obtain for La 
KB (x x) Гале, — me SL 45 — x) 
and for Lp: 
К) (x — 2’) = – ie {ГГ лье, — πιο (iT E + Р) τὴ Ag A (x — x) 


и’и“» 

With 
A (x — x’) = — A(x’ — x) 
9A(x—x) AA (x – х) 


Ox θα’ 


the supplementary condition for Ф gives the following identity 
Күк — x )zAKT (a —x) А 


where А is defined by help of formula (4.1). 
After substitution of (5.12-13) in (5.15) we get 


κως KI КӘ ко) 


for arbitrary isotensors in К and KO), 


(5.14) 


(5.15) 


(5.16) 


It may be easily seen from (5.16) that the Lagrangian. Ly is the true Lagrangian of our 
quantized equation, if from the two possibilities in (5.3) and according to (5.11) the anti- 


commutators are chosen. 


For commutators we obtain instead of (5.15) 


K(x — х) = — AK? (x' — x)A 


(5.17) 


and the Lagrangian £ describes boson fields, satisfying our generalized Dirac equation’. 
The relation (5.13) fully determines the commutation relations between a and 


(ge Let us introduce 
Ф) — РФ 
where, assuming 
(poo) = pu) 
PPD — РР) = 0 


(5.18) 


(5.19a) 
(5.19b) 


* The problem of *spinor-bosons" satisfying the Dirac equation will be investigated elsewhere, 
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we may put for e=—1 


T h ‚ ΑΠΡ τ 


үр № TE 
The subsidiary condition implies 
Ф = АФС) = АФ (5.21) 
Thus from both 8-spinors (5.18), we shall take into account only the 8-spinor DP, 
Using 
Ф — (perg (5.22) 
we obtain 
{ФСР), Qo Ην ΚΤ (5.23) 
and because 
[P= К®] 0 (5.24) 
from (5.195) we have 
{DP OF} =0 (5.25) 
Similarly 
(OZ, OE} = 0 (5.26) 


Using (5.22) it is easy to show that 
(DR, OY} = (РОК ӘР) z (5.27). 
or, by virtue of (5.19а) and (5.24) 


BP, ФУ?) = DER) (5.28) 
Calculating (5.28), а 
(1.6. Бут. РӘ) = „үг? F2. PH (5.29) 
we get finally 
OH), oC?) = {VDE FICE De] (Are, me a A(x — 22) 


where є,, is a unit isovector paralell to #,,. 


6. The Jauch Field 


The field which obeys the commutation relations (0.5) was first proposed by Jauch 
in 1954. 
We introduce the transformation 


D> ФВ = ВФ (6.1) 
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which transforms (0.5) in the conventional anticommutation relations for the Dirac field, 


namely 


WR), valx’)} = 0 (6.2) 
Let us insert the Majorana representation of the y,-matrices (Majorana 1937), 

for which 

у= Ср рї (6.3) 


We have 


1 Y = LPE ET 
@ = — Ф = — (y, — ys) (6.4) 
ү? p yg - v 


The commutation relations (0.5) can be written by help of (6.4) 
Ky (х л less (їз — 5190) (1,9, T; — т)} лв A(x — х) (6.5) 


and similarly (6.2) 


Κάδία -α)-- (rf ur — m))4g* A(x — х) (6.6) 
Putting in (6.1) 
В = exp (- als (6.7) 
we obtain from (6.6) 
BK P(x — x')B = exp (us) - K Q(x — η (6.8) 
Because | 
exp (т, Г.В) + τα = я IEEE (6.9) 
if 
cosh В = im (6.10) 
] — p? 
after the scale transformation 
Φ-» ) Ф 
Vine (6.11) 


we obtain from the Dirac commutation relations (6.2) the Jauch commutation relations (0.5), 
where the transformation (5.7) is an isorotation in the (1,4)-isoplane. 


The transformation (6.7) do not change the isovector #.,, chosen for the Dirac and 
Jauch fields by means of (0.1). 
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r ^ . $ 
The full 3-parameter group, which leaves the mass isovector components (0.1) inva- 
‘lant, is 


y = ap + ἓψ |а|®— |02 = 1 (6.12а) 
or, using 8-spinor notation, 
; un i 
Ф: =.exp meo А x Grat Ф (6.12b) 
It may be divided into two subgroups Q and J: 
li 
О: Ф = ехр (9:7 Ф (6.13) 
J: RO — exp {— 4 (0171 ΞΡ от) Г5}Ф (6.14) 


The l-parameter group Q leaves the free antisymmetrized Schwinger Lagrangian 
(Schwinger 1948) invariant: 
dEn = ivy, Ἔ πι)ψ 5 Ly (y n ar m)y* (6.15a) 
or, using 8-spinors as in Sect. 4, 
Lo = aal 8,13 + т)Ф (6.15Ъ) 


О can be used as an electric charge gauge-transformation group, not changing the Lagrangian 
and the conventional commutation relations (6.2). 

The 2-parameter group J leads to a generalized Jauch field, described by the commuta- 
tion relations (0.5), (0.5b). Ш we transform (5.6) by means of (6.14), we obtain after a scale 
transformation 


Dm Φ 6.16 
| — of -- 02 f 

that in (0.5b) 
0 = 01 10 (6.17) 


The Jauch field differs from the Dirac field only by another choice of the commutation 
relations. Thus we see that in quantum theory a 4-spinor field is not determined by the 


equation uniquely. 


7. The Second Isovector ε΄, 
Ву (4.15), (4.8) and (5.6-7) we have (for f =0) 
d KE DU ЕЕ Dy + (ГЭьт,е „+ m)? (2.1) 


MUS 
where 


eme (7.1a) 
m 


5 
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The Lagrangian (6.15) of Schwinger is obtained from (7.1) by putting (0.1) and 


Е=0 D=(0,0,1) (7.2) 


It is easily seen that (6.15) transforms under (6.7) to the following Lagrangian of Jauchl 
field: 


Lo = Φίτα + Г; oral Zut + m)® (7.3)! 
Using (4.8) апа (4.11) we have 


Е = (0,0,0) ἢ-- (0,0,1) (7.3a)) 

From (7.1) it may be seen that the Lagrangian of our generalized field is determined] 

by the choice of e, and Εν = (E, D), or, taking into account (4.8), by the choice of 6, 

and FR... | 

Let us consider the connexion between €,. €, and Е. Using (3.22а), (3.20а)) 
and (3.21), we obtain 


к, = — Ве „Е, (7.4)) 
and 
gus. s εδ. (7.5) | 
where 
C = — P (7.5a) | 


From antisymmetry of E,,, we obtain the relation 
А “A / e А = 
65,6, = C-6 Fé, = 0 (7.6) 


uv zu 


: : uere A : eu т С / Р i d 

which by virtue of (7.1) is equivalent to (0.2). Therefore, є, and €, are orthogonal. As it 
“ зи 

is proved in Appendix I that 


в„Ё „= 0, (7.7) 


р | 
the isovector c. 
HO 


is orthogonal to the isoplane described by the isotensor Bs 


w v 
It is necessary to mention that if the invariant E. Е vanish, Е. „апа E,,, describe 
34V 7 SUP sur jur 
two orthogonal isoplanes. Because for each isoplane in 4-dimensional space, there is only 
, = : κ H 
one orthogonal isoplane, we see from (3.16), (3.21) and (7.7) that є and є, are parallel 
3 
to the isoplane described by E,,„,. Taking into account the identity (7.6), we conclude that 


. , , . . . H 
the isovectors é,, and e, (or х; 6) describe a 2-dimensional orthogonal coordinate 
system on this isoplane. 


Comparing (6.15) and ( 


4.8), we obtain for 6 in Schwinger's case the relation (0.3), 
and in Jauch's ease — (0.4) 


The isoplane, orthogonal to c. 

H 5 H A 

charge isoplane. In Schwinger’s с 
isoplane. The rot 


and €, (described by means of ЕЁ.) is the the electric 
ase it is, in accordance with (0.1) and (0.3), the (1.2)- 


ations in the electric charge isoplane describe, in accordance with (6.13), 
the electric charge gauge-group О. | 
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9. Conclusions 


_. The theory is determined by a pair of orthogonal isovectors in Minkowski’s isospacez 
: | Ee : і 
Because (see (3.8), (7.1a)) always e, > 0, we have the following two invariants with respect 
to the C'-rotation classes of the Lagrangians: 


И (8.1) 


euclidean electric charge isoplane 
and 


є ==! -ᾱ- eg ЕЙ (8.2) 


(pseudoeuclidean electric charge isoplane). 
The Dirac-Schwinger choice (0.1)—(0.3) belongs to (8.1). 

It may be shown that we pass from type (8.1) to type (8.2) after, e.g., by the following 
interchange of the Lagrangian (6.15b): 


Lo = ФГ (Г.Э, ль + m)® (8.3) 
which leads to the commutation relations 
| 


πα wey} — 0 
{p,(x), ya) = — ὑδαρία — x')+ A(x — х) (8.4) 


Both cases (8.1) and (8.2) admit the introduction of the following generalized commuta- 


tion relations 
(φαί). wll = αδαρία — 2}4{α — х) 
(φαί). gll = В5ав(х — м) (ах — х) (8.5) 
where o, Û are complex and 
ques p (8.58) 


The form (8.5) of the most general commutation relations is written for the choice (0.1) 


of the mass-isovector components. 


9. Final Remarks 


In this paper the quantum theory of 4-spinor fields, invariant with respect to C x С” 
‘is developed. From our point of view, it may be interesting to discuss the following two 


problems: 
1) The connection of the 4-spinor theory with the 8-spinor theory, which describes, 


as we suppose, the free baryon field. It may be mentioned that the supplementary condition 
Ф=АФ (9.1) 


can be treated as а baryonic number neutrality condition, 


2) The introduction of the interaction Lagrangians, invariant with respect to € x С” 

We intend to investigate further these problems elsewhere. 

The author would like to thank Professor Rzewuski for his advices and many, helpfu 
discussions. 


APPENDIX I 


. . , . . . . . . 
If a certain isovector €,, exists, which satisfies (3.21), then the mixed invariant. 
5 


obtained from the two isotensors En. E,,, vanishes identically. This theorem may b 
simply proved by obtaining first the result that the isotensor 


Ομ = ЕщЕ, (A.D) 
is diagonal. 
Indeed, the relation (A.1) can be written 


1 


Giu T 9i T M OR (А.2у 


Let us assume u 52 v. Thus, if оу > o, and Ш # 01 #1, U #0, v, the isotensor component 
(А.2) is equal to the following sum of four terms? 


С» = = { 


7 f ; 
9] LE ме: бубен E seus St E uo Gran δ μα] Li 9i VE wo Crowe, x 


X Ее, + E uo E voe E suo.) (А.З) 
Using the antisymmetry property of Е u» and єє, we obtain the result that the first! 


term in (A.3) is equal to the second, and similarly, the third to the fourth. We have 


and 


ο. SCH b. d du [ῃ (A.5) | 


he form (A.5) is not invariant under C’-isorotations. Thus 


tu =t = FSP G. (A.6) 
Because 
Sp Cie ei D eue (A.7) 
we have | 
Guy — ZE (A.8) 
From (A.8) we obtain 


Gi, Fa Gy (A.9) 


5 We do not use here the Einstein summation convention, 
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ind, using (3.21) 


ο EG ms EE =0 (А.10) 


a 5HO 20 


On the other hand, it follows from (A.8) that 


SIE ἔπε, (A.11) 


and, if such а v exists for which €,, 70, we obtain (3.16). Now we show the validity of (7.7). 
Using (7.5) ава (A.9), we get 


= Si RT, : s D = / 
EEL, = С eL ELE, = С° 6,0, = O° ga 


ШААН) 


/ id / x 1 
CO сбив, = С Eye Е 


4. 
220 


(А.12) 


and if relation (3.21) is valid, the right side of (А.12) vanishes and relation (7.7) is proved. 
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MEASUREMENT OF THE MEAN LIFE OF THE FIRST 
EXCITED STATE OF Na? 


Ву B. Амввойу, A. Fauprowicz, A. Jasinski, J. ΚΟΜΝΑΟΚΙ, Н. LANCMAN, 
J. LUDZIEJEWSKI 


Polish Academy of Sciences, Institute of Nuclear Research, Warsaw 
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The mean life of the first excited state on Na?? was measured using resonance fluorescence 
technique. Ne?? obtained from Nei? (n, у) Ne?’ reaction by irradiation neon in the reactor was 


used аз a source. The value of the mean life was found to be т = (1.5 ze x 1072 sec. 


1. Introduction 


In view of the possibility of obtaining information on f-v angular correlations by 
using nuclear fluorescence technique [1], [2], we have undertaken the investigation of the 
decay of Ne?3, which is a suitable isotope for this purpose, as suggested by Lewis and Cur- 
tis [2]. There is a marked discrepancy between the results of measurements of the mean 
life of the first excited state of Na23 published by various authors: Booth её al. т= (1.5 + 
+0.3) x 10-10 sec [3], Krone et al. 1—10-12 — 10-13 sec [4], Rasmussen et al. т=(1.8 + 
+0-4) x 10-1? sec [5]. Therefore we have undertaken to measure this mean life once more 
applying the self-absorbtion method [6], [7]. 

Results concerning the decay scheme of Nei? will be published separately. 


2. Experimental 


'The experimental arrangement is shown in fig. 1. №23 was obtained from the Nei? 

(n, y) №23 reaction by irradiation of spectrally pure Ne in the reactor “EWA”. The thermal 
n . 

neutron flux around the activation chamber was 1077 = The gas was circulated by 
a centrifugal pump driven by an electric motor. The voltage on the motor was stabilized. 
Since the intensity was dependent on the pumping speed, care was taken to work in an interval 
where this dependence is approximately flat. Therefore the intensity of the Ne?? source 
was constant within 1%. The brass source chamber had a cylindrical shape and dimensions: 
gi 30 x 40 mm and wall thickness 1 mm. The possible source inhomogeneity did not exceed 
5% along the axis. The source intensity during the measurements was about 100 uC. 
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А 2" x 2" NaJ(Tl) crystal mounted on EMI 6097B phototube was used. The pulses 
| were analysed by a 100 channel pulse height analyser. The half width of the spectrometer 
was 12% for 661 keV line of Cs137, 

As the source intensity was small a thick scatterer had to be used. The sodium scatterer 
and absorber were prepared in thin wall plexiglass containers by melting the sodium in vacuum. 


The thickness of the absorber was 4 cm. The dimensions of the scatterers were as follows: 
Ф..= 15.5 cm, d 7.5 сш, thickness 4.0 cm. The comparison scatterer was made 
of crushed MSIE του the comparison absorber of aluminium plates. The dimen- 
sions of the comparison scatterer and absorber were the same, as those of the sodium. In 
order to match both scatterers for Compton scattering and pair production, a small amount 
of paraffin was added to the magnesium scatterer. The matching of the scatterers was checked 
with Au!’ and Pr? sources. The absorbers were checked by comparison of the attenuation 
of the intensity of the Compton scattered y-rays from the magnesium scatterer. The relative 
sensitivity of the NaJ(TI) crystal for different points of the scatterer was measured using Au!%, 


3. Results 


The measurements were made in the neighbourhood of the reactor. Sometimes, in 
spite of thick asphalt-lead shields, some fluctuations of the background, in intensity and 
shape of energy spectrum, were observed during the measurements. It was established that 


this fluctuations with good approximation could be considered as linear vs energy in the 
resonance region. Taking this into account the spectra obtained with the sodium scatterer 
were normalized to this obtained with the magnesium scatterer. The corresponding curves 
obtained in this way are shown in fig. 2. The rather large background in the 440 keV region 
is caused by the high energy part of the y-radiation near the reactor. 

Fig. 3 shows the difference curves obtained by subtraction of the spectrum measured 
with the magnesium scatterer and aluminium absorbers. The fluctuations of experimental 
points at the left edge of each spectrum are predominantly statistical in origin (partially 
they would be caused by random errors in the placement of the scatterers). In the normaliza- 
tion procedure, mentioned above, both spectra have been cut off along a straight line connect- 
ing their right edge and the valley between the compton edge and the photo-peak. This 
is why the spectra are asymmetric and they go to zero at the left edge. The shapes of the 
spectra marked out by the experimental points are consistent with the shape of the 440 keV 


line of Ne?3. 
5 
The measured intensity ratio (see appendix) is — = 0.51 + 0.04. Using the equation 2 of 
/ E 


the appendix, with 7,,— 305° [8], [9], /,—3; h=% [10], ΠΠ, [12], [13] we obtain 
К — 0.27 + 0.045. From this the mean life Sg R "to ВЮ 
This is in disagreement with the result obtained by Booth οἱ al. [3] in a similar experiment 
as ours and in good agreement with that obtained recently by Rasmussen et al. [5], who used 
the 440 Key y-radiation from Na?? (p, p^). From their results т = (1.8 703) x 10-1? and the 
known partial lifetime for electric quadrupole transition (E2) —5.6 x 10—19 sec. [14], [15], 
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amplitude ratio д =0.056 and from angular distribution of resonantly 


scattered y-rays they found ὃ to be rather positive. 

The mean life calculated from the unified model [17] and from the single — particle Weiss- 
kopf formula (for M1 transition) is 4 x 10-13 sec, which is not too far from the experimental 
value. 

The authors would like to thank Dr Z. Wilhelmi, chief of Nuclear Physics Laboratory, 
for the kind interest shown in the present measurements, Mgr Eng. J. Aleksandrowicz and 
all members of reactor crew for many facilities in performing of the measurements. 


APPENDIX 


Taking into account the geometry shown in fig. 4 it is seen that the number of counts 
of resonantly scattered (assuming isotropy) y-rays is given by. 


χ 
не { 
ο) ο КУЛЕ «qm 
Fig. 4 
S(D.K) = cons ATAT en N 2 
S (D.K) = const Д R? €(5. 0) е F(D + д, К) (1) 
where 


dt, = dxdyda 
ат, = одоабаф 


(— Ка)" 


F CH = Nr LE 
(К)? = т! (m+ 1)% 


1 F(q, К) was tabulated by S. Ofer and A. Schwarzschild [7]: 
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Τε is the effective temperature of the sodium scatterer and absorber. 
n is the number of sodium nuclei per cm? 
ε(ξ, 0) is the relative sensitivity of the NaJ(T1) crystal for various parts at the scatterer 
€ (С, о) = 2 (6, о) e-wte—ov. є' (C, о) was measured experimentally when the scatterer was removed. 
μ is the total nonresonance absorption coefficient of the scattering material. 

If we put in formula (1) instead of 4=9(х, y. a, C, о) and R= R(x, у, a, ©, о) the values 
q=q(x, y, a, C, o), R =R(x, у, a, C, ϱ) averaged over x and y, which is a good approximation 


for our geometry, we obtain by integrating over Ф 


| «aca, = € (6, o) F(D+ д, K) e-"45, (C, о) 


6 SDE 


5 d 
EET atm 700 = 0) 


(2) 


= 


ER 
Θα 15 the cross section surface (for a given fixed a) of Ще source chamber seen from the point 


5 x c 
(С, 0) of the scatterer. The curve — vs К numerically calculated from equation (2) for our 


So 
geometry is shown in fig. 5. For comparison — = F(D, К) for the case of a thin scatterer is 
: So 
also plotted (the lower curve). 
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ÜBER DIE QUANTENTHEORIE DER ANTIFERRIMAGNETIKA 
Von LEON KOWALEWSKI 
Institut fiir Theoretische Physik der Adam Mickiewicz Universitat, Poznan 
(Eingegangen am 16 Februar 1961) 


Es wird mit Hilfe der Methode von Holstein-Primakoff die Dispersionsformel fiir Fer- 
rite mit einer inversen Spinellstruktur berechnet. Wegen der magnetischen Überstruktur erhält 
man eine anisotrope Dispersionsformel. Im Falle gleich grosser Spins bei den Oktaeder-Gitter- 
plätzen ist für Antiferrimagnetika E ~K?. Es wird gezeigt, dass die Ergebnisse der korrigierten 
Arbeiten von Wonsowski und Seidow und von Kondorski, Pachomow und Siklós sich gleich- 
wertig sind und mit den Ergebnissen der vorliegenden Arbeit übereinstimmen. 


I. Einleitung 


In der letzten Zeit erschien eine Reihe von Theoretischen Arbeiten (z.B.:Wonsowski 
u. Seidow 1954; Kondorski und seine Mitarbeiter 1957; Kaplan 1958) über antiferrimagne- 
tische Körper mit einer Spinellstruktur. Das Wesentliche dieser Struktur ist ein aus Sauerstoff- 
Ionen 0— gebildetes, kubisc!-flachenzentriertes Gitter, in das die Metall-Ionen Me**, Ее? 
eingebettet sind, wobei sie Tetraeder-Zwischenräume (A), beziehungsweise Oktaeder-Zwi- 
schenrüume ( B) einnehmen. Es wird vorausgesetzt, dass die dreiwertigen Eisen-Ionen in den 
Zwischenriume (A) einen Gesamtspin o, die dreiwertigen Eisen-Ionen in den Zwischen- 
räumen (В) den Spin S und die zweiwertigen Metall-Ionen in den Zwischenrüumen (В) den 
Spin s besitzen. 

| Die Rolle der Sauerstoff-Ionen in der Wechselwirkung zwischen den magnetischen Ionen 
berücksichtigen wir, indem wir Superaustauschintegrale einführen. 

Die Ergebnisse der obenerwühnten Arbeiten widersprechen sich gegenseitig. Die von 
Wonsowski u. Seidow (1954) abgeleitete Dispersionsformel für Spinwellen besteht aus zwei 
Energiezweigen (Ё) die beide linear von der Wellenzahl К abhängen, wenn man nur die 
Wechselwirkung zwischen den (A) und (B) lonen berücksichtigt (die 4-4 und B-B 
Wechselwirkungen sind verschwindend klein im Vergleich mit der А-В Wechselwirkung). 
Die Dispersionsformel aus der Arbeit von Kondorski und seinen Mitarbeitern (1957) besteht 
gleichfalls aus zwei Zweigen, die beide vom Quadrat der Wellenzahl abhängen. Sowohl Won- 
sowski u. Seidow als auch. Kondorski mit seinen Mitarbeitern stützen ihre Rechnungen auf 
der Theorie des polaren Metallmodells. Kaplan, der in seinen. Rechnungen die von ihm 
erweiterte Methode von Anderson (1952) benutzt, nimmt an, dass im allgemeinen die Spin- 
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wellenamplituden bei den sechs Arten von Gitterplätzen im Spinell ungleich sind. Bei dieser 
Annahme erhielt er eine aus sechs Energiezweigen bestehende Dispersionsformel, von denen 
viere vom Quadrat der Wellenzahl K abhängen und zwei gleich grosse konstante Grössen sind. 

Wenn man gleiche Spinwellenamplituden bei den Gitterplätzen A und gleich grosse 
Amplituden, wenn auch nicht dieselben, bei den Gitterplätzen B voraussetzt, da erhält man 
nur zwei Energiezweige. 

Kaplan behauptet, dass die Annahme gleicher Amplituden in den Bewegungsgleichungen 
zur Formel EK?, und dieselbe Voraussetzung im Hamiltonian, im Falle gleich grosser Spins, 
zur Formel E-K führt (siehe Kowalewski 1961). Kaplan sieht den Fehler von Wonsowski 
u. Seidow darin, dass Sie im Hamiltonian unrechtmässig gleich grosse Amplituden annehmen. 

In der vorliegenden Arbeit wird mit Hilfe der Methode von Holstein u. Primakoff die 
Dispersionsformel für Spinwellen berechnet. Es wird gezeigt, dass, im Falle gleicher Spins, 
die Energie der Spinwellen auch bei der von Wonsowski u. Seidow eingeführten Annahme 
von gleichen Amplituden vom Quadrat der Wellenzahl K abhängt. 


II. Der Hamiltonian und seine Eigenwerte 


Die Berechnung der Erregungen nach der Theorie des polaren Metall-modells (Bogo- 
liubow 1949) ist völlig gleichwertig der Methode von Holstein u. Primakoff (1940), wenn 
ausschliesslich Spinerregungen vorliegen. Der Hamiltonian von Bogoliubow (1949, Formel 
4.109) geht in diesem Falle über in die Diraesche Formel für beliebige Spins. 

Der Hamiltonian hat also folgende Form: 


ο AMA Dn" 5, UI 


J( Ji» №) Superaustauschintegral zwischen den Gitterplätzen / und fs; 


2» 


ζω 


τ : : dë Ae Et 

5; Spinoperator beim Gitterplatz f in = 

f Index der „magnetischen“ Gitterplátze. 

Wir gehen zu den Bose-Operatoren b über: 
97 = 5, (1), + bre(n;)] 
3 WT CN EE ea eee 
57 = ing V Sy [g(n,)b, — bj (n, )], 
ὃς (9L*L с 


A n DNK ч : 
ev К и E 5, ny P bj by. 5, i. Ki SS 5, D тын 


np = l für f welches einen Oktaeder-Gitterplatz 1 oder m bezeichnet; 
ue = — 1 für f welches einen Tetraeder-Gitterplatz r bezeichnet. 

Das Spinell-Gitter ist aus sechs Translationsteilgittern zusammengesetzt. Die Vektoren 
der elementaren Translationen haben folgende Form (siehe Abb. 1): 


> (9 жем Ke ES а = а 
وة ,)6+8( >= ,)6.45( = ت‎ = + 8), 


e 
€,» €y, €, — Einheitsvektoren, 


Abb. 1. Elementarzelle eines Spinells. 


Lagevektoren der Tetraeder-Gitterplätze: 


Кп) p re 2s 0, Or, i» 0, 
3 
> > a E 557 
Ry NE пау, Or, 4 елге δε, 
j=1 


М, МА = 


N Gesamtzahl der ,,magnetischen” Gitterplätze. 
Lagevektoren der Oktaeder-Gitterplätze: 


R(n) = R, + 6), δι = 


a 


01, E 8 (36, zi 5e, gir Зе.), 


у; 


= > > a > 
p (n) РЕ? IU F m: Om, = 8 (3e, al 7еу 


> 


A η (ee). 
Wegen der Bedingung der Quasisättigung erhält man aus (1): 


ЕВ + Σ [20 Σ de vi — 28 Y Ilr, ἢ — 2s D Ar, В, + Σ 5»; TEOS 
ТА 1 т 
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E Ze Y, (L7 τ 25. m)] 192 [25 Σου m’) — 20 У) Im, r) +255, Din 


— 2e У! Ir, т); by, ρα 25 Σὲ οι ON a 


(r, т’) (m, m’) 
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-. 2/65 > Jr. Db E b b) — 9V'os > VS m)(b,b,, + b,5,) — 
(r, Ὁ (г, т) 
— 2/55 У) Kl, m)bjb, — 2/55 У) Kl, π) δ) (3) 
(1, т) (l, m) 
wo 


Ey Co $0? У) οτι, г) – 1.52 У) ДЬ, 1) — 455 Σ πι. т.) ep al? 


LE d lys my nts 


+ as p J(r, m) — sS >) J(L m). 


1, т 


-— OE TA 
Wenn wir im Hamiltonian (3) Лу, = - ЛЬ, fa) setzen, ist er bei der voraussetzung 
o = 5 = s identisch mit dem Hamiltonian von Wonsowski (W.S.1.7.). Bei der Annahme 


о = $ = s = 1 ist er identisch mit dem Hamiltonian von Kondorski (K.10) oder (K.24). 


Um den Hamiltonian (3) zu diagonalisieren, zerlegen wir die Operatoren b in eine 
Fourier-Reihe. 


Vektoren des inversen Gitters: 


3 
SSC 1 * l 5 
їд 7 ] 
= =e NIS У, 
j=1 
> gy 
misc de 


Nehmen wir an, dass die Fourier-Transformaten der г Operatoren 6, und b, gleich gross 


und eleich b sind. Nehmen wir ebenso an, dass die Fourier-Transformatender Operatoren 
5 a P 


Ge b, 4 δ, Om, gleich gross und gleich ©, sind. 
Die ο e haben dann folgende Form: 


(4) 


Die Operatoren b,, c; erfüllen folgende Beziehungen: 


^ 


Δ Wi ^ Ax Ak ^ 
bib, — b, b, = сус, — € C, = Oy. 


Den Hamiltonian (3) kónnen wir jetzt in der Form 


^ 


Н = Е,4 2 АБД Σ Bb e USE bees D Dee, (5) 
a A 
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chreiben, wo 


A, = 2e У) ът) [1 — exp iK,- (В, — R,)] - 25 У) J(r.l) — 2s У) J(r. m), 
r l m 


= Ὁ 
σ = mI > => 
ae -V =) [νο A. An, Dep i (Ё, = R) + 
p=1 l 


iE ys 2 (τρ. m) exp iK, . (R,, — Rn) | 
|5 p ЛО, D) [1 — exp i£, (В, — RJ] + 


SC 3 Im, т) [1 — exp iK; (Rm Й„)]+ 5 ἃ Im, ἢ) | » 
m 1 


= у: ехр К. (В, — Bnp | +s E Л, т) [ πα ехр ТА kK — δρ == 
—g Σ J(l,r) --σ Y J(m, т) (6) 


Wir vernachlässigen den Einfluss der Kristalloberfläche. 

Den Hamiltonian (5) diagonalisieren wir mit Hilfe der Transformationen: 
TE os P 1 
— (дар), (у == = 


v2 ү? 


О, = (1+ w)-* (way, д? da); Pici (= ру) + Рэд), 


b; = (Ree νη. 


R= (Lu) (Gut wda)  Sa= (1 — w) (рід — WPa): 
Den Parameter w bestimmen wir aus den Gleichungen (Tessman 1952): 


Wir erhalten die Eigenwerte: 


Е — const — 3 Y [A,+D, VA; D) – 483] + Enna t 2) Baue: (1) 
A A A 


В = А.И (А, + D)? — 483], 
pod e Dot V(A; + р)? AR), 
Ny, = 0, l, 2, ..., пал = 0, 1, 2. ... 


(Та) 
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Kondorski mit seinen Mitarbeitern erhielten die Grössen Ag, Bk, C,, D, (siehe: 
Formel K. 18). 

Wenn wir die Ausdrücke (K. 18) korrigieren, indem wir die Spinellstruktur berück- 
sichtigen, und in den Formeln (6) о = $ =з = 1 setzen, da erhalten wir folgende 


Beziehungen: 


A, Аһ  B,-y2B, С = 75 Въ Dp = D. 


Man sieht dann leicht, dass die Dispersionsformeln von Kondorski (K. 17) mit den 
Formeln (Та) (für о = S = = 1) identisch sind. 

Wenn wir in der Arbeit von Wonsowski und Seidow in den Transformationen (W. S. 1.8) 
die Normierungsfaktoren korrigieren, in 4, Bi, D; (У. S. 1.9) die Rechnungsfehler 
beseitigen, und die Superaustauschintegrale Лл, durch die Integrale — 4 J (fi, fo) ersetzen, 
so wird die Dispersionsformel (W. S. 1.17) identisch mit dem Ausdruck (7), wenn 
στ S = s ist. 

Mann sieht also, dass die korrigierten Arbeiten von Wonsowski und Seidow und von 
Kondorski, Pachomow und Siklós zu denselben Resultaten führen (für Spins vom Werte 1). 

Wenn wir uns beim Rechnen zu den nüchsten Nachbarn und zu langen Spinwellen 
(kleines K) beschrünken, so erhalten wir eine folgende Dispersionsformel: 


Ej, = QU + Qa? K? + Da Kr)? + ΩΘα:Κ«Κ- (8) 
wo 


(5 Я 
о@ e (34) 1 7 SE 728 Blon ES ПЛ] +6/4в(5 + 5 —o) = 2JBBÛ, 


1 
1 e 
ον = + =. (УЈвв — 4oJ44) + (96%) —1[330(/$ + Vs) Лав + (уЛвв --4oJAA)Y(B — 30) ], 


1 
πο (=) 


32 B 
8 ο. M un ze Јвв 
248 == (5 = 5) + Увв KS (3%) 150 Лав + (30 = В)Јвв = D (За) 
16 Jan 
a = Vo? + S? + 52 255 — 40/5520. 
Be 540-2155, ος ο ο ο κος 
ΡΕ Superaustauschintegrale zwischen nähsten Nachbarn 4 — A, A — B, 
B — B. 


Die Dispersionsformel (8) mit den Koeffizienten (8a) ist anisotrop. 


Die Anisotropie der Dispersionformel (8) ist durch die magnetische Uberstruktur 
verursacht, 
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Im Falle S = s (6 42S) erhalten wir: 
ФО эж Бу = 5) J a5; ΄ 
1 1 
Q® = o) ep Qe zm 
) 


1 


( ыы 
QP = + 3 (5/вв oJaa) + 57 (08 E EE o ER (68) 
Die Formel (8) mit den Koeffizienten (8b) ist isotrop. 


Wenn o=S=s ist, so erhalten wir: 


Ἡ T 
Ore ου RE 0, 


Т 
©(Јвв + 2JA4). 


33 
в 
а SJAB T 


ste irr 
Ө? = zg 4B = 194a + 2J n). (8c) 


In der Arbeit von Wonsowski und Seidow kann man dieselbe Abhängigkeit E ~k? 
erhalten, wenn man in ihr: 1) die schon Oben erwühnten Fehler beseitigt; 2) wenn man auf 
richtige Weise die Kristallstruktur berücksichtigt und die Rechnungsfehler in der Náherungs- 
berechnung der Dispersionsformel beseitigt. 


Wenn о = 25 = 2s ist, d.h. im Falle eines Antiferromagnetikums, und wenn man nur 
die A—B Wechselwirkung berücksichtigt, so erhalten wir die lineare Dispersionsformel: 


Ед = ὅρη = ‚ "ok, (9) 


Herrn Prof. Dr S. Szezeniowski danke ich für wertvolle und nützliche Diskussionen. 
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КРИТИЧЕСКИЕ МАГНИТНЫЕ ПОЛЯ ТОНКИХ 
СВИНЦОВЫХ ПЛЕНОК 


Ежи Шимашек 
Криогенная Лаборатория ПАН, Вроцлав 
(Поступила в редакцию 16 февраля 1961 г.) 


Были измерены критические магнитные поля для тонких плёнок свинца 
в температурном интервале, равном около 1,2 К ниже критической температуры 
свинца. Измерялись плёнки толщиной 3,15 . 1079 см — 1,16 . 1074 см. Для тонких 
плёнок фазовые переходы из сверхпроводящего состояния в нормальное или 
наоборот, в присутствии магниного поля — обратимы. Для тольстых плёнок пере- 
ходы обратимы лишь вблизи критической температуры Тк, тогда как для температур, 
более отдалённых от Тк, переходы сопровождаются появлением гистерезиса. 
Уже при температурах близких Тк получено максимальное значение гистерезиса. 


$1 Введение 


Свойства сверхпроводников, которые по размерам сравнимы с глубиной про- 
никновения магнитного поля в массивный сверхпроводник (0—10 9 см), отлича- 
ются от свойств массивных сверхпроводников. 

Впервые исследования тонких плёнок сверхпроводников были проведены неза- 
висимо A. И. Шальниковым [1] [2] и Эпплярдом, Bpucroy, Г. Лондоном и Мизе- 
нером [3] [4] [5]. 

Шальников, в результате исследований тонких плёнок олова, свинца и таллия 
толщиной 5. 10-7 см— 3: 10-5 см, пришёл к выводу, чт овеличины критических 
полей для тонких плёнок значительно больше, чем для случая массивного металла. 
Измерения для свинца проводились лишь при температурах ниже 4,2° К. 

Проводя подобные исследования для тонких плёнок ртути, Эпплярд, и со- 
авторы нашли значительное увеличение критических магнитных полей по срав- 
нению с массивным металлом; кроме этого, в своей работе они приводят зависи- 
мость отношения критического магнитного поля тонкой плёнки Нк к критическому 
магнитному полю массивного металла Нкм OT толщины плёнки d. 

Эта зависимость имеет следующий вид: 


(1) 
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b — равно приблизительно 1,2: 10 5 см при низких температурах и зависит OT 
температуры так же, как глубина проникновения магнитного поля. 

Зависимость (1), подтверждённая также с помощью опытов Эндрю [6] и Лок- 

ка [7], справедлива лишь для толстых плёнок, причём величина b не сохраняетс 
при изменении сверхпроводника. 
Нк 
Нкм 
более тонких плёнок, пользуясь феноменологической теорией сверхпроводимост 
Ф. Лондона и Г. Лондона [9]. 

Гинзбургом и Ландау [10] была разработана другая форма макроскопическо 
теории сверхпроводимости, в которой, кроме свойств массивных сверхпроводни- 
ков, была разработана проблема сверхпроводимости образцов малых размеров. 

Из теории Гинзбурга-Ландау следуют особенные свойства тонких сверхпро- 
водящих плёнок, находящихся в магнитном поле. 

Согласно этой теории, плёнки, в зависимости от толщины, можно разделить 
на два класса: dd и а> dy, где dg является некоторой критической TOJI- 
щиной. 

Для тонких плёнок (d < dy) переход из сверхпроводящего состояния B нор- 
мальное или наоборот, в присутствии магнитного поля, является фазовым перехо- 
дом второго рода. Для тонких плёнок справедлива следующая зависимость: 


Неудачной оказалась попытка Гинзбурга [8] найти также зависимость дл 


Ho. wt e 
Ho η 9) 


Для случая толстых плёнок (d > ἐκ) переход из сверхпроводящего состояния 
в нормальное или обратно является фазовым переходом первого рода, причём | 
для плёнок с 4>4к выполняется зависимость: 


Нк 8 
Hot V (3) 


Заварицкий, проводя измерения на тонких плёнках олова [11], таллия и индия: 
[12], полностью подтвердил выводы, следующие из теории Гинзбурга-Ландау. | 

В настоящей работе приводятся результаты измерений магнитных критических ` 
полей для случая тонких плёнок свинца. 


ὃ 2. Экспериментальная часть 


С целью получения температур выше 4,2° К, был построен криостат, пред- 
ставлённый на рис. 1. Похожий криостат был предложен Заварицким [13]. При- 
менение автоматического фотоэлектрического усилителя для питания нагревателя 
позволяло поддерживать температуру с точностью 0,0005 °К в течение довольно 
долгого промежутка времени. 

Для измерения температуры был использован угольный термометр сопротив- 
ления, сделанный из радиотехнического сопротивления типа Аллена-Брэдлея. 


555 


Рис. 1. Криостат для температур выше 4,2 AR 


1 — сосуд Дюара с жидким гелием, 2 — перевернутый дюаровский сосуд, 3 — Экран из медной 

фольги, 4 — нагреватель, 5 — медная полоска, соединяющая экран с жидким гелием, 6 — вен- 

тиль, Ти 9 — угольные термометры, 8 — катушка из свинцовой проволоки. На рисунке нет 
внешнего дюаровского сосуда с жидким азотом. 
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Сопротивление использовавшегося термометра, изменявшееся в температурном! 
интервале 6,2°K—7,2°K на 9900 омов, измерялось с помощью высокоомного по-. 
тенциометра. Градуировочная кривая термометра контролировалась перед каждым. 
экспериментом в двух точках: температуре кипения жидкого гелия под атмосфер-. 
ным давлением и точке перехода массивного свинца из сверхпроводящего в нор-. 
мальное состояние (7,22°К). Для исследований употреблялся свинец высокой. 
чистоты, не меньше, чем 99,999%, что было проверено спектральными и химичес- 
кими методами. 

Образцы были получены путём конденсации свинца в вакууме на плоских, 
оптически полированных стеклянных подложках. Четыре платиновых провода, 
пропаянные сквозь подложки, использовались в качестве токовых и потенциаль- 
ных проводов для измерения сопротивления плёнок. 

Толщина плёнок определялась весом испарённого металла. 

Для определения перехода образцов из сверхпроводящего в нормальное со- 
стояние измерялось изменение сопротивления образца в зависимости от темпера- 
туры (рис. 2), или величины внешнего магнитного поля. 


De є не ar > - 

Рис. 2. Изменение сопротивления образцов разных толщин в зависимости от температуры 
К — значение сопротивления образца 
KR, — значение сопротивления образца в нормальном состоянии 
Г — температура в °K. 


В первом случае определялась температура перехода Tk плёнки, во втором — 
значение критического магнитного поля Н к При определённой постоянной темпера- 
туре 7. В качестве Ty и Нк принимались значения, которым на определённых кри- 
вых перехода соответствовало сопротивление, равное половине сопротивления 


09% 


бразца в нормальном состоянии (вблизи температуры перехода). Сопротивление 
бразцов измерялось с помощью низкоомного потенциометра. 

Магнитное поле получалось с помощью катушек Гельмгольца (Н >>250 эрст), 
или электромагнита. 

Магнитное поле было установлено всегда параллельно плоскости образца. 


$3 Результаты 


Было исследовано поведение в магнитном поле шестнадцати тонких свин- 
цовых плёнок толщиной 3,15.10-6 см — 1,16.10-* см. Исследования проводились 
в температурном интервале AT = Тк — T, равном для некоторых образцов около 
1,2°К (Тк — температура перехода плёнки в отсутстви магнитного поля, Т — тем- 
пература измерения). 


| 7000 
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3,63-102 
413-106 
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100 6.49 10-6 


667108 
993109 
111-105 
14410° 
462-1075 
219-1075 


323.105 
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Рис. 3. Зависимост критических полей Нк от температуры AT. Пунктирная линия — значения 
критических полей Нкм для массивного свинца. 
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Температурный интервал, в котором сопротивление образцов менялось от 
‘значения К, (в нормальном состоянии) до нуля (в сверхпроводящем состоянии), 
у большинства плёнок не превышал 0,03°K. | 

На рис. 3 показана зависимость критических магнитных полей плёнок Нк от 
температуры ДТ. Kak видно из рис. 3, для плёнок толщиной меньше, чем 1,62. 105cm, 
сохраняется характер зависимости Нк = f (AT) во всём температурном интервале, 
причём критическое магнитное поле пропорционально АТ”. 


0830: 


250 


Рис. 4. Кривые перехода из сверхпроводящего B нормальное состояние в магнитном поле для 
плёнки толщиной 3,23 - 1075 см 


Н — напряжение магнитного поля в эрстедах 

Цифры над кривыми перехода обозначают величины AT — Тк-Т 
О — результаты получённые при увеличении поля 

€ — результаты получённые при понижании поля 


Область зачерченная — область гистерезиса. 


Для плёнок толщиной больше чем 1,62. 10-5 см на кривых зависимости H, = 
/ (AT) наблюдается излом; влево от излома характер зависимости Hg = f (AT) cox- 
раняется таким же, как и у плёнок тоньше чем 1.62: 10-5см; напротив, вправо от 
излома зависимость Hp = f (AT) похожа на зависимость. характерную для массив- 
ного металла (пунктирная линия на рис. 3). 

Было найдено, что для плёнок тоньше 1.62: 10 


-CM фазовые переходы из cBepx- 
проводящего состояния в нормальное, 


происходящие в присутствии магнитного 
поля, — обратимы. Наоборот, для плёнок более толстых, чем 1,62. 10-5см, фазовые 
переходы обратимы только вблизи критической температуры плёнки; 
с определённого значения 


начиная 
ЛТ, наблюдается явление гистерезиса: значения магнит- 


ных критических полей при переходах из сверхпроводящего состояния B нор- 
мальное больше, чем значения при обратных переходах (рис. 4 и 5). Как видно из 
рис. 4, область гистерезиса увеличивается с ростом AT, crpe 
значению. В таком поведении состоит разница 
и оловом [11] и индием [12] с другой. 


МЯСЬ К постоянному 
между свинцом с одной стороны 


Рис. 5. Кривые перехода для плёнки толщиной 7,47 · 10 ? см 
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Рис. 6. Кривые перехода для плёнки толщиной 9,83 · 10-9 см 
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Было замечено, что для тонких плёнок переходы из сверхпроводящего в нор- 
мальное состояние и наоборот являются крутыми (рис. 6). Для более толстых 
‘пленок крутые переходы характерны лишь вблизи критической температуры; 
начиная с определённой температуры ЛТ, сопротивление растёт круто от нуля до 
некоторого значения меньшего, чем значение сопротивления образца в нормальном 
состоянии. При дальнейшем увеличании магнитного поля сопротивление очень 
медленно растёт, вплоть до величины сопротивления R,, а затем не претерпевает 
последущих изменений. 

Takoe поведение более толстых плёнок может быть связано, как нам кажется, 
с возможностью возникновения промежуточного состояния для свинцовых плёнок. 
Более определенно этот вопрос могут решить дальнейшие исследования. 

Автор выражает глубокую благодарность А. И. Шальникову за постоянную 
помощь во время выполнения этой работы в криогенной лаборатории Московского 
Государственного Университета. 
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ТЕМПЕРАТУРНАЯ ЗАВИСИМОСТЬ ГЛУБИНЫ 
ПРОНИКНОВЕНИЯ МАГНИТНОГО ПОЛЯ 
В СВЕРХПРОВОДЯЩИИ СВИНЕЦ 


Ежи Шимашек 
Криогенная Лаборатория ПАН, Вроцлав 
(Поступила в редакцию 16 февраля 1961 г.) 


На основании измерений критических магнитных полей тонких плёнок свинца 
приводится температурная зависимость глубины проникновения магнитного поля 
в сверхпроводящий свинец. Экстраполируя получённые значения к T=0°K, автор 
получил значение 09—(4,5 + 0,2): 10-8 см. 


Ilo cux пор в литературе встречается мало работ, посвящённых измерению 
глубины проникновения магнитного поля в сверхпроводящий свинец. 

На основании опытов по измерению сопротивления тонких свинцовых проволок 
Понтиус [1] определил величину глубины проникновения магнитного поля с точ- 
ностью до порядка (10-5 —10-6см). 

Стиль [2], исследуя магнитную проицаемость сверхпроводящих свинцовых 
шариков радиусом порядка 10-3см и используя зависимости, вытекающие из те- 
ории Ф. Лондона и Г. Лондона [3], определил величну глубины проникновения 
магнитного поля при температуре 4,22°K. Согласно его данным она равна (1,3 4 
+ 0,3) - 10-*cm. 

Локк [4], используя теорию Ф. Лондона и Г. Лондона для полученных им 
экспериментально магнитных моментов тонких плёнок при температуре 4,24°К, 
вычислил глубины проникновения для пяти свинцовых образцов. 

Принимая температурную зависимость глубины проникновения магнитного 


поля, которая следует из теории Гортера и Казимира [5]. 


б 
à = ewm (1) 
T 
SG 
Тк 
где д, — значение глубины проникновения при T 2 OR 


T — температура 
Тк — критическая температура 
последний получил для Du величину, равную (3,9 = 0,3) - 10-9cw. 
(565) 
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Симон [6], используя результаты измерений поверхностного импеданса на 
частоте 9200 мггц и формулы из теории Ф. Лондона и Г. Лондона, представил 
графически температурные зависимости глубины проникновения. 

Глубины проникновения магнитного поля, получённые из измерений активз 
ного и реактивного сопротивления, отличаются между собой во всём температур- 


1 ab 04 42 01 005 


ТЯ 


Температурная зависимость глубины проникновения магнитного поля для свинца. . Henpe- 
рывная линия — результаты автора Пунктирная линия — построёна по формуле (1) 


ном интервале. SE значения глубины проникновения для T — 0°К 
даёт значение Ó, = 5,4. 10—6см. 


На основании измерений магнитных критических полей тонких плёнок свинца 
3 
выполнённых автором [7], также можно вычислить и дать графическое представ- 


ление её изменения с температурой. Для этой цели удобно использовать зависи- 
мость: 


a τ. 2۷6 а Θ) 


вытекающую из теории Гинзбурга и Ландау = 
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где: Нк— критическое магнитное поле тонкой плёнки 
Нкм-— критическое магнитное поле массивного металла 
д — глубина проникновения магнитного поля 
d — толщина плёнки 
Необходимо добавить, что зависимость (2) справедлива лишь для таких тонких 
плёнок, для которых фазовые переходы из сверхпроводящего в нормальное со- 
стояние и наоборот являются переходами второго рода. 

С помощью непрерывной линии на рисунке представлена температурная за- 
висимость 0(AT), получённая автором (AT=T,—T). Пунктирная линия представ- 
ляет зависимость (ЛТ), вычислённую из формулы (1). Экстраполируя получённые 
результаты измерений для д к температуре T'—0^K с помощью зависимости (1), 
автор определил глубину проникновения магнитного поля в сверхпроводящий 
свинец. Полученная величина равна: A 


д = (4,5 + 0,2) - 10cm. 


и находится между значениями, получёнными Локком и Симоном. 
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ERR CRITICAL PROBLEMS FOR MULTILAYER SLAB SYSTEMS 
By А. KUuszELL 
Polish Academy of Sciences, Institute of Nuclear Research 


( Reecived March 1, 1961) 


In the work the solution of the critical problems for slab multilayers systems are given 
by the procedure developed by K. M. Case for one-velocity Boltzmann equation with isotropic 
scattering of neutrons. In all considered cases the solution of Boltzmann equation is reduced 
to the solution of one dimensional Fredholm type integral equation with an additional critical 
condition. 


Introduction 


In the report the critical problems for slab multilayer systems are solved by the procedure 
derived by K. M. Case [1]. The following types of problems are considered: a) multiplying 
slab with infinite reflector, b) three-zone core with infinite reflector and c) slab with finite 
reflector. In all cases symmetry with respect to the central slab is assumed. In all cases the 
solution of the Boltzmann equation is reduced to the solution of the Fredholm integral 
equation with an additional critical condition. The solution of the Fredholm equation cannot 
be expressed in an analytical form, an approximate method using the spherical harmonics 
treatment of the kernel is proposed. 


1. The solution of the Boltzmann equation and the properties of eigenfunctions 


One velocity Boltzmann equation in plane geometry in the case of isotropic scattering 


has the form 
1 


H E i LZ = if Wa u) ам’ (1.1) 

—1 
where index i denotes respective region, и is the cosine of the angle between neutron velocity 
and x-axis, J; is the neutron mean free path and с; is the number of secondaries per collision. 
Following K. M. Case’s [1] results we can express the solution of Eq. (1.1) in the form 


x 


x 
v», 
H 


1 x 
P(x, u) = oe Wi Φίμ, v;) + бе d (ui, — v;) + { Aj(v)e ' Diu. v) dv (1.2) 
—1 


(567) 
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where 


y 


ga T λιμ)δίμ — A (1.3) 


y — 


Φ(μ. v) = a 


In the later considerations all integrals involving terms of the type will be un- 


derstood in the sense of the principal value. The function Я; (и) is given by the formula ı 
Alu) = 1 — cy tan A-1 y (1.4) } 
and +», are the roots of a transcendental equation 
i 
cy; tan h-1— = 1 (1.5) i 
Vi 


It is well known that for all c; v; not belonging to (— 1, 1) and Ф, (м, + »;) reduces to ı 
the form 
Vi 


Фи, x») = 2 (1.6) 


ZU FU 


The functions 4; (v) and the constants а; and 6; can be evaluated by means of the bounda- 
ry conditions depending on the problem considered. 
As it has been proved by K. M. Case the functions of the form (1.2) 


are complete in the 
space of solutions Eq. (1.1) 


D 


It is easily seen that the functions Ф, (и, v) satisfy the following "orthogonality" relations 
1 
C —C& m 
utu, v) Фи, аи = nAiv)AkGó(v — 5) + ist (1.7) 
1 
С; — Ср 
äu +») Pla, v)du = vOv, + ν)) (1.8) 
Wi 
2S 
Cj — Cy jv 1.9 
1 EV doizk E 
РЕ 
КУ E vi).(u. + vy) du = 
Se aL T Sy ‚)-2 | ΝΙ(-- νὴ for i=k | 
| z 0) E j 00 | ( 
and 
| μμ. + η) Dy (и, F νχ)άμ = + Mr MER (1.10) 
H 2 Vi ns Wi 
-1 


IR. H : А 3 : x 3 
Putting i=% we obtain the orthogonality relation derived by K. M. Case for the function 
defined in the same region. 


Using Eqs. ((1.7)-(1.10)) and the Poincare-Bertrand theorem we obt 


ain the following 
formulas necessary for furthe 


r analysis 
1 


1 1 
een [omas Фа 2. D 000%, v)dy 9 
= Ἐπ e 


1 
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al 1 
D fin, v)du f QGp du. η)άη = v [ * 2 seo? 00) + 
—1 


=! 


1 
o co | ndn ` ‚ 68 — 6 ОР 
ae 5 of aor vNa(v)Q(v) + - ne SES ( (1.12) 
| 21 
The formulas (1.11) and (1.12) are valid for arbitrary functions Q(v) which satisfy the 
ólder conditions. 


2. The critical problem of a slab with an infinite ref lector 


Let us consider an infinite slab of multiplying medium (c>1) of the thickness 20 with 
reflector of infinite thickness on both sides of the slab. The Cartesian coordinate system 
ith the origin in the centre plane of the slab and x-axis perpendicular to this plane is chosen. 

The solutions of the Boltzmann equations will be denoted as follows: for the slab 
var, и) and for the reflector a(x, u) for х>о and W(x, u) for х< —o. The boundary 
onditions in our case are as follows: 


a) The solution must be continuous at the boundary between the slab and the reflector 


V (о, u) = Pole, и) 
Y,(-0 м) = #3(—0, и) (2.1) 


‚ b) The solutions in the reflector must vanish at infinity 


lim 45 (х, u) = lim Vx, u) = 0 (2.2) 


x>00 x>—00 


The Boltzmann equation in plane geometry Eq. (1.1) is invariant with respect to the 
transformation of (x, u) into (—x, — u). The symmetry of our problem implies that the 
boundary conditions are also invariant with respect to this transformation. Therefore the 


solutions are also invariant 1.е. 
V (s, -- μ) = Ῥη(--α, — p) (2.3) 
and 
P(x, u) = P3(— x, — p) (2.4) 


From the above conditions we can conclude that it is sufficient to consider only the 
region x20. The extension to the region «<0 is straightforward. Using the above results 
the boundary conditions can be expressed in the ultimate form 


7(0, u) = 7,0, — p) (2.9) 


Ψη(ο, u) = Palo, и) (2.6) 


lim Ps(x, u) = 0 (2.7) 
x-00 
From the conditions (2.5)- (2.7) the coefficients of the eigenfunction expansion can: 
be evaluated. Inserting the function (x, и) given by Eq. (1.2) into the condition (2.7) we: 


obtain that 


Zz»)s0 - forc 9-0 (2.8) 
and 
"m (29) 


Multiplying both sides of Eq. (2.5) by иФ, (и, v) and иФ\ү(и, т.) respectively, andi 
integrating with respect to м over the interval (—1, 1) we have the following relations:: 


1 1 Р 
f Dru. v) V0. u)du = | UD, (и, v) Y (0, — u)du (2.10) | 
-1 -1 
and 
: ` 
J “Фи. ev). Aide = | ифи. + о), — дар (2.11) 
-1 —1 


As the functions Ø, (и, v) and Ð (и, +») are complete in the interval ( — 1. 1) the relations : 
(2.10) and (2.11) are equivalent to Eq. (2.5). 
Using the orthogonality properties of the functions Ф\(и. v) and Φι(μ. +r) we obtain 
from (2.10) and (2.11) the following relations: 
A,(v) = Aq(— v) (2.12) 
and 
a, = b (2.13) 


Further relations can be deduced by multiplying Eq. (2.6) by the functions иФь(и, v) 
and uG,(u. +v), respectively, and integrating over u €(—1. 1): | 


\ 
6 : 
ῬΝ(ν)ε ^" Α.(ν) =| uPs(u, v) (о. udu (2.14) 
-1 
Sek i 
Νου "" a, = 1 Balu, νι) (0. u)du (2.15) 
-1 
SIE. 1 
№, (— т) бе ^" = l ИФ. (м, = νι) Vo, и)аи (2.16) 


-1 

Nsy(v) and Nal +5) are defined by Eqs. (1.12) and (1.9). 
Eq. (2.14) for »—0 and Eq. (2.15) determine the coefficients of the eigenfunction 
expansion in the reflector region by means of the solution in the multiplying region. Eq. 
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.14) for v<0 and Eq. (2.16) using Eqs. (2.8) and (2.9) can be written in the form 


di 
| ban, — fe, даи = 0 fo ν»0 (2.17) 
—1 
nd 
1 
Tun, — 9 (o, ди = 0 (2.18) 
== 


These equations with the conditions (2.12) and (2.13) can be reduced to the integral 
quation for the function A,(v) defined in the interval (0,1) and to the additional condition. 
he constant a4 cannot be evaluated in the problem considered as the system of Eqs. (2.17), 
2.12), (2.13) and (2.18) is homogeneous with respect to 4,(v) αι and Ьу. Thus all the consid- 
red equations can be divided by the constant αι. Ultimately the integral equation has the 
orm 


Е n | _ 20 
т 7 Ei Co ) А | GE 4 ; 
Nalu) Y(u) + 5 IE ραν Υ(ν)ά) 
0 


C1 == Со 


"m е 
un SE Sa | (2.19) 


1 


where а new unknown function given by the expression 


[^] 
e A, (vje 
үр) = 200°" (2.30) 
y 
is introduced. 
The additional condition 
_ 201 zx. d 
[АД [АЛ li 
АРТ Б, (2.21) 
БР И Wat Vo Ур 7» 


is to be considered as a critical condition. The solution of the Boltzmann equation for the 
reflected slab has been reduced to the solution of Eq. (2.19) with the critical condition 
(2.21). These equations give the angular distribution of neutrons in the multiplying medium 
as proportional to an arbitrary constant. The angular distribution is directly obtained by 
means of the above relation. 
Let us consider the limiting case when &>0. From the definitions we have 
a) lim Λο(ν) = 1 
cz 
b) lim νο = оо (2.22) 
DE 


c) lim N) = 40) 


ει 
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Putting these values into Eq. (2.19) we obtain the integral equation of the form 


1 
| _ 20 
lv 
г с у Де ο 
Alu) Τιίμ) 4f e se "IT Y(v)dı 
0 


Же Tg И | 
= — posu. —w)e ™ + @,(u, поля] (2.28 


The critical condition (2.21) tends to the identity 0—0. | 

The equation (2.23) is identical with the equation derived Ьу В. Zelazny [2] in the case 
of a bare slab. 

The equation (2.19) is a singular integral equation which can be reduced to the Fredholn 
equation by the method given by N.I. Muskhelishvili [3]. This Fredholm equation haa 
a form (see App. A) 


1 1. 
> ur Ίο πα У S(T 
Y(u) - I | ve BY γί) Е τ᾽ | un А Ala 


2T(u) , Mä ш (n + v) "cru 
^ 
£ aa | ПИФа) W (n) 294 
= доди) — эз | PAVE a, (2.24 
0 
where 
^ л? ә 
AY(u) Au) + ам | 
S(u) = s; ; EN 5 (2.25} 
Luna: Tac! T | : 3 πμ] 
T(u) = Xis(u) + T CCo? + a E = aiul (2.26) 
2 Re | 
Dee? л? нео 2.97) 
Wu) = Xielu) MA As) + I" Sr ER (2.27) 
and 
Су — С, =. £ 
Фа(и) = — 3 E : kou, νι) οι + D (u, — ve ση (2.28); 
1 


The function Au) is defined by Eq. (А. 9). As the index of the singular integral equation: 
(2.19) x —0, the equations (2.19) and (2.24) are equivalent without any additional condition. 

The solution of the Fredholm equation (2.24) is not known in an analytical form. One: 
of the possible approximate methods of solving this equation based on the spherical harmon- 
ics representation is given in App. B. 
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3. The critical problem of a three-zone core with an infinite reflector 


Let us consider a system consisting of a central slab of the thickness 201, characterized 


| ‘material parameters сү, /, surrounded by two symmetrical slabs of the thickness 05—01 
#scribed by material parameters с», lą. The whole system is immersed in an infinite re- 
Ῥοῖοτ, the values of material parameters of which are denoted by с; and Js. 
| The origin of x-axis in the Cartesian coordinate system lies in the plane of symmetry 
the system. 

After a similar discussion as that of the Chapter П, the boundary conditions can be 
pressed in the form: 


V (0, и) = V,(0, — u) (3.1) 
Ῥηίοι, ш) = Vp, и) (3.2) 
(05, ш) = Vos. и) (3.3) 
: г 
lim Miz, u) = 0 (3.4) 


Where the indices 1,2,3 denote the successive regions, beginning with the Central 


Ёр. 


Let us write the general solutions (1.2) of the Boltzmann equation in a more convenient 
orm 
X—n—1 ~~ On 
P(x, 0) = aye "^" Φίμ,ν))}-ὄ.ε"'" Фи, —») + 
1 
χ-θῃ-ι D dmi 
+ | E ^^ x) B(—»e ” EH Ф (и. v)dv (9.5) 
-1 
fore πια 
nd - 
_х=@ х—оз 
ψη(α, м) = азе I" Dalu, νο) + Бе "* Фу, — >) 
1 E 
+ f {AH + DN ais >” Dalu dv (3.6) 
—1 
here 


Ἱ {πε О 37 
το) = for v<0 η 


is the conditions (3.1) and (3.4) are identical with the similar condition in the Chapter 
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II the solution can be readily obtained in the form 


d; 
Ai») =e ” B,(v) ER 
A ' 
а =е "b, (3 
where 
d, es = (3.1 
1 
and 
B,(v) = 0 (3.11 
5, =0 (3.12 


Multiplying both sides of Eqs (3.2) and (3.3) by w®,(u, v) and u®,(u, +25), respectivel 


and integrating over u ei 1, l) we obtain the system of equations: 


1 
M 
ОН -»e'x(-»[-— [ Au. ν) Pilo, u)du (5.13 
—1 
1 
_ 4 | 
Nato) Ayo " X(v) + В, (— v) {= bese: V)P3(0., u)du (3.14 
-1 
E 
No(s) as = J Osu, və) (01, μ)άμ (3.155 
sl | 
_ ἀν ' 
№(— νε "b, = | иФ,(и, — νε) Y^ (o, , u)du (3.16) 
σα 
-d ` 
Na(va)e "ra, = MPU, vy) 4(0, и) (3.17% 
—1 І 
1 
N δν = | иФ,(и, — v) Vos. μ)άμ (3.18) 
-1 
where 
Cy О 
d, = SR (3.19). 


Eliminating the expansion coefficients with the index 2 from these equations we cani 
reduce them to the form: | 
а, 1 


1 
e * f p®,(u, fe, аи = | nui. v) (os. и)ди (3.20) 


—1 dc 
for 0 
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1 
iB 
иФ,(и, — ν)ψη(οι. du =e » | иФ,(и, — ν}9' 
à Jos ο (3.21) 
fo v>0 
ET 1 
e^ f ифи, sde, даи = | иФ,(и, +v) Paler udu (3.22) 
-1 -1 
ith the “orthogonality” relations of Chapter I we have. 
1 1 
q (КО 9) (os, μλάμ = sell As) + 2— 2» | A) и а» 9 Фу, κ) 
2 ER с : 
4 $ [ 8 
| (3.23) 
ing Eq. (3.23), the Eq. (3.20) has а form 
1 d τ 
Eu [uds „1 e = 
В+ [ по due Ze | иаи ан — as B Φιν, v) 
= 2 : 
(3.24) 
his equation is solved in App. A. The solution has a form 
1 
| Eu 
: ^ e v = 
413») = Кз | ‘as | иФь(и, v) (01, μ)άμ — аз = 7 2 (v, νο) (3.25) 
| 3 


=I 


here the operator Gë is defined by Eq. (А. 25). Similarly as in the Eq. (3.23) we have 


1 1 
BL A 
Laun = gll (05, μ)άμ =— Ra [559 du (3.26) 
- 0 
for ν--0 


hich with the Eqs (3.11) ава (3.25) gives 
1 
| иФ,(и, — v) Pag, u)du 


; 


d. 
сасы -- d . dg BG 
о, [E к, m ηΦη(η. u) (от ап + аз 7 2 Ф. (и, v3) 
3 


0 τα (3.27) 


Jenoting the term in the function V (o,, м) connected with the continuous spectrum by 
P (u) the function Y (о, u) can be written as follows 


di 
la: u) = V, (u) x au. νι) + bye Ἢ du, = νι) (9.28) 
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Using Eqs (3.9) and (3.28) we have Eq. (3.27) in the form 


1 
ELA — yy V. (u)du + 
-1 
3 1 
15 
Joan di ν | 
ает | EE Ra τ. 0 фыр дату = fo) + ago) G2: 
a u 
| St 
where 
. = DL 
А) = 3—-3 E ^ QGi(—»,»,-- Φι(-- >, — νι) + 
SS 
qp ceno μις \ 
Ee ee) E Kaye в — E ” Ø (u, νι) + Ou, — n) | | (3.9 
2 1-8 = | | 
1 
P (сз — с)? _® мам o, - SI 
gi») = 2с» pe J ut ; Fas {P3(u. v3) (3.33 
0 
Similarly, Eq. (3.22) can be expressed by 
1 1 
+ B eS) sty ο тс. 
е dg νο) V (оу, u)du — 2—2 ^ B didum: 3 ЕЛУ ›„)4и (3.33 
E e Vy Hs Pos 
Inserting Eq. (3.25) into the last equation we obtain 
1 
*[ Cg — Co Vo 
| ^ I usu, + νο V^, (o, μ)άμ = ———3 —23 
EYE, 2l, 9) A. μ)άμ OB νη 
-1 
1 1 | 
BS 
+ an Ф,(и. +v) K. | x nD (n, uy) (о. та Ca 68 с 
с, Ш 540, E Vg A K (lau. ШР (0.2) а) — аз ο Ф. (и, уз) Аи (3.38 
—1 —1 


which can be written also as 
1 
= d: 


d, 
dot Laux + νο) ε(μ)άμ + 
1 


1 


Ps Фи. + vy) Rs SH пФ. м) (т) аи = a + b+ y + аз (3.34 


vhere 
d, 2d; 
EF EC па uir 
2 Vor Vy νο + | 
1 
en "E Tu, = SL А С νε " + Q(u, — ZI (3.35) 
0 
1 
aaa” Е 
dek = FRE Vs τ THE. bn, + va) Къз (Pau, v3); du (3.36) 
б : 
With the conditions (3.8) Eq. (3.34) can be reduced to 
& bı + уаз = Г. [Βι(ν)] (3.37) 


where the functionals Г, [B,(v)] are given by 


: 2d 
T2 с — ο -— 
Γι[βι(ν)]}--ᾱ.” : Jan ij » D + νο) + 4,(— v, + a” ἄν + 


(ο 


1 
г Г.С 
+A | Ra | sn в) ΠΠ δη) + 


0 
+ Φιίη,-- E dvd» du (3.38) 


Dividing both sides of Eq. (3.37) by аз and denoting 
Eu 


iS (3.39) 
3 
b 
ge E (3.40) 
we can rewrite it in the form | 
1 Е 
Mët Se? TRAO) — ys} (3.41) 


Comparing the above equations for quantities with plus and minus indices, respectively, 
we obtain the condition i 


1 1 A 
a (30) y RETI {Г_[Ү,0)] – у-} (3.42) 
Pob XE 
The similar procedure applied to the equation (3.23) gives it in the form 
1 
Ww jl H 
№0) Tall + 25 af рҮ) = du / Au, v) Υπ(μ)άμ = = — Lao - Σὲ gp) 
ge id Hi 
0 б 


(3.43) 


where 
1 1 2d, 
QU, رھ‎ BS | BE ا _ ا‎ 
ή 0 
d : _% . ' а, 
πα ο ο ENTIS | e "Фи, & + Φιίη, | аат ἀμ 
2 μον μι у, - 


0 0 —1 
| (3.44) 
It can be easily seen that the kernel (и, v) is а Fredholm type kernel. Therefore Eq. (3.43) 
can be reduced to the Fredholm equation by treating it with the operator defined by (A.25). 

As the index of the singular integral equation (3.43) is equal to zero, the reduction to 
the Fredholm equation does not involve any additional condition. Therefore the solution 
of the considered problem is equivalent to the solution of the Fredholm equation 


1 
^ ^- 1 γα, 
Ττ(ν) + Koy [ου ») T E kV Е x Ze yo | (3.45) 
0 


with the additional condition (3.42) which plays the role of a critical condition. 
The approximate method proposed in App. B. can be used in obtaining the solution 
of this equation. 


4. The critical problem of a slab with an finite ref lector 


Let us consider a system consisting of a central slab of the thickness 20,, characterized 
by material parameters су and i, surrounded by two symmetrical reflectors of the thickness 
02—01 described by material parameters с» and l. The origin of x-axis in the Cartesian 
coordinate system lies in the central plane of the slab. 

We can solve this problem in two ways. The first one is based on the previous consider- 
ations. Putting c,—0 in the problem considered in Chapter Ш we obtain the identical 
system. However, the passing to this limiting case is not trivial because the operator Kos 
changes its class so that a new form of the additional conditions is generated. In the second 
method the system with finite reflector is treated independly from the beginning. In the 
Chapter II it was showed that both methods of treatment of finite systems are equivalent. 
In order to avoid tedious algebra in further analysis the second way is chosen. 


Using the symmetry properties of considered system the boundary conditions can be 
written as follows 


1(0, u) = V (0, — u) (4.1) 


τοι. и) = Ῥηίοι, и) l (4.2) 
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па 
Ἔρος») =й for ue (4.3) 


The solution of Boltzmann equation in our case will be taken in the form (3.5). 
The condition (4.1) wbich is identical with the condition (3.1) can be reduced to 


d; 
Al =e ”Вү(›) (4.4) 
nd 
di 
αι-ε ^b, (4.5) 


Let us write the condition (4.3) in the explicit form 


E & 2 1 
азе "Ф,(— M, v3) baal и, = va) H | A,(rje "Φιί-- u, »dr + | BA u, —v)dv 
0 0 


O (4.6) 


Using the properties of the function Ф„(и, v) we can write Eq. (4.6) in the form 


1 d; 


1 
c, [ vB. c, [ve "As(v Eo 

Λο) Balu) + ES | v = ya d ys A ! dv — азе ^ ,(u, — vy) — byby(u, νο) 
6 


(4.7) 


The relation (4.7) has a form of the equation discussed in details by K.M. Case [1] and 


R. Zelazny [2]. 


The solution of this equation has the form: 


7 ба р \ 
^ v A (v t 
Ρεμ) = — Ks E E y+ " dy Базе "Фи, — vy) + θοΦοίμ. vs) (4.8) 
0 
where the operator К. can be written аз [2] 
а λο(μ)/ (u) 
Ka{f(u)} = ry 
[49 (ш) + E " 
1 - 
1 с пХ Qf n)dn (4.9) 


у 2 . 
Xz (м D = "Те | Lou = Gol [η = μ] 


0 
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where 
1 1 ` d 
XzQ) ER уан | тг (4.10) 
ET 
0 
and 
λο(η) + iz Ca 


E 
λο(η) — i э °27 


The singular integral equation (4.7) has the index х=1 so that additional condition is 


generated: 
1 e 1 d, 
ЛЕ --- | 9 y de E 
HH PL ае ** D,(n, — νο) + b, (n, νο) + 3 AL Їу|ар = 0 
2 ER 
λο(η) — i5 CoN : 


0 


The condition (4.12) can be rewritten in the form 


34 IURI E 
e з ха + Db, + = E " h(v) 4 (ν)ν = 0 
ὃ 
where 
1 ^ 
ي‎ (x ۰ EA 
t 
en) 
К: 
1 - 
go (Xeno дал 
λο(η) — 5 ομη 
` 2 
and 
H - 
h(v) = nAz (тап 
Е - Fein [η -- >] 
0 


- From the condition (4.13) we have 


1 
b = e Ж do — ο 
2 p 2 9 В 


0 


а 


Е η (v) Ag(v) dv 


(4.14) 


(4.19) 


(4.16); 


(4.17) | 
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Inserting Eq. (4.17) into Eq. (4.8) we obtain 


1 
$ C5 aA T / 
Ви) = = DE № " A,(v) E Dp T dAn. 5] 2 = 


d 


I3 Lo (оын, --ν) — 5 Pall ») (4.18) 


Using the formulas (4.17) and (4.18) we can express the function V. (от, и) by 4ο(ν) 
and a, only. j 


1 


a ads ка: 
NOM CD) e. ^ e "+ 5 ve улу) Φ,(μ. — νο) + 
SH: 
1 1 
„© Κη α \ 
+ [ лано. ν)άν — f " Ф.(и, — v) E s Clou, — v) — B Φ.(ν. el = 
0 0 | 
Er 
Ba -2 l h(n) 
ck zelt, n A(n) EZ ОЙТ. D,(v, „| an} | dy (4.19) 
0 
The last expression may be written as follows 
H 1 1 
Vos. н) = as) + | Add, vide" [ Обь v), 0v (4.20) 
0 0 
where 
1 
= аа = а EU R 
Ф(и) = yu. vy) — В enu (0200) Bu jn E "un, — v) К oso es 
0 
mo 1 (4.21) 
Be te | 
and 
1 1 
" we 
| Qi. элд» = 55 | Фи, — »)h (ve? Αιλάν — 
0 0 


ЕЯ 


- а = а. 
micas), EI bt, 
2. 2 e A 2 fe " As(n) | l = 2 Φ.(ν, „| dn) D (u, — v)dv (4.22) 
) 0 


It can be easily seen that the kernel Q(x, y) is of the Fredholm type. Multiplying both 
sides of Eq. (4.2) by the functions u®,(u, v) and μφιίμ. +»,) respectively, and integrating 
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over ш € (—1, 1) we obtain the system of equations 


dı 1 


νΝμ(ν)ε ° A,(») = [ UP, (и. v) Polo. u)du (4.23) 
-ᾱ 
for у> 0 
1 
— УМ 1 (>) В: (у) = | uD (u. --ν) 'Po.u)dyu (4.24) 
SH 
for у> 0 
1 
Νι(νι)αι = | uu. vy) Vo. u)du (4.25) 
-1 
4 E 
— Ne "b = | uD (u, — n) (0. u) du (4.26) 
Using the above relations we can derive from Eqs (4.4) and (4.5) the following equations 
1 2d: а 

| nus 2) (0, du = -е f аФи. — FG u)du (4.27) 

Б for » > 0 

and 
: 06. 
] ud. (и, νι) Paloi κ) = — е" | иФ\у(и, — vı) (01. u)du (4.28) 
= —1 


With Eq. (4.21) the last two equations can be written as follows 


1 1 1 


^ 


а» | μμ. »)P(u)du + | vun) “Ὁ Qv, νι]άν + | As(r) | uD (u, νι) Q(v, u)dudv 
А . Cy ^ 
0 


0 e 
24 К - 
=—е ^» fas | UP, (и, — ALIMEN > | vd (v) 2 — 2 Φι(ν. — νι) ἂν + 
-i 0 2 
en 
dÉ КО | HO (и, = v) Об, u) dudy (4.29) 
D -1 
and 
1 1 1 1 1 


| uD (и, »)Dlu)du - EZ (u, v) db Jann, ndy 4 d Aaly) / ИФ (и, Or, u)dudv 
0 


2d, d 


2а, 1 = - * 
AN d J utin. — Du) du + {| иФү(и, — ») | An) Dau, т) + 
—1 —1 0 
1 1 
+ | 4:0) f ntu. it, юл (4.30) 
0 —1 


for v>0 
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Dividing the both sides of Eqs (4.29) and (4.30) by a, and using the “orthogonality” 
properties of the function Ф,(и, v) [=1, 2 we obtain Eqs (4.29) and (4.30) in the form 


к : С; — 6 126 
- Το(ν) Es ν(Φι(ν,νι) * Φι(ν, — 4)) + 
0 


2d; | 
+ | о. и) (ous. nie  Φιίμ, — J in| dy 
=й 


T 


fa 24; 
= | иФ(и) (ai pre O J du (4.31) 
= 
and 
1 1 d 


5 - Ίο τν fi Ve Oa 3 = > 
v Na (+) ο(ν) 4 Co п gis lié ο) ἀμ’ Co = Eu [655ε SC 


0 


1 1 
ή 2d; 


a | Ya | 001.0 | E о = J dudy 
0 -ᾱ ў x 


1 
2d; 


= — ED | ο δα = ‚| Яша (4.32) 


е. 


where 


(4.33) 


Eq. (4.32) can be rewritten in the form 


1 


1 
В — с n > 
кызуу) a Tan mz, 
0 


e 


iy 
2d, 
2» _ | иФ(и) ο. ντε * d, — d a Ss 
= | 


where 
2d; 


2d, 


ο — б ue 1. ] --- 
use) = = 5 = vm 5 | Ξ | nQ(u. 1) |90» +e 2 Out, »| dyn (4.35) 


This is the singular integral equation the dominant part of which is identical with the 


equation considered in App. A. Using the operator К, we can reduce Eq. (4.35) to the 
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Fredholm equation of the form: 


1 
У. (+) + Let Ω(μ. v) Υγίμ)άμὶ 


21 | ) 
0 


1 


5 3 2d, 
== E | иФ(и) | уте 7 Φι(μ. > 2 (4.36) 
0 
without any additional condition. 
So the equation (4.36) together with the critical condition (4.31) uniquely determines 
the solution of Boltzmann equation. 
It is clear that method of approximate solution proposed in the App. B as well as other 


approximate methods can be used in this case, too. 


5. Conclusions 


In the report the following critical problems been considered: 

a) multiplying slab with infinite reflector 

b) three-zone core with infinite reflector 

c) slab with finite reflector. 

The solutions of all this problems have been reduced to the solution of the one Fred- 
holm integral equation with the additional critical condition. The solutions of respective 
Fredholm equations are not known in analytical form. | 

However, the approximate numerical procedures can be used. Such а procedure involv- 
ing the spherical harmonics representation has been proposed. 

The proposed spherical harmonies method is different from the classical one in the 
point that only the transport part of the solution is expanded into spherical harmonics. 
It is worth to note that in all considered problems the same form of the operators 


^ 


К» and the functions Ху (и) is met. Therefore the operators and the functions can be 
tabulated independently what the problem is considered. 

The generalization of the method for the systems with the greater number of layers 
on for the non-symmetrical systems should not lead to the basic difficulties. The increase 
of the number of layers gives only the increasing of the number of resulting Fredholm 
equations. 

For instance, the solution for the three-zone reactor with finite reflectors can be 
reduced to the solution of the system of two Fredholm equation with one additional condition. 

In all considerations it was facility assumed. the none of the сүз is equal to the unity. 
However, when one of су tends to the unity, the discrete eigenfunctions tend to the form 
given e.g. by Van Kampen [4] all results are still valid. 
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APPENDIX А 


Reduction of the singular integral equation to the Fredholm type 


The general theory of singular integral equations given by N. I. Muskhelishvili [3] can 
e used in solving the equation of the type 


1 Фо) Y(») e 


Ο(μ)Υίῳ – = | и) (A.1) 
0 
where 
Ош) = Аи) Aut) + F аси (A.2) 
and 
Plu) = ли (АЗ) 


According to Muskhelishvili's method we introduce а sectionally holomorphic function. 


Mrs [ SORUN (A.4) 


det ZI, po 


Using Plenely formulas we obtain from Eq. (A. 1). 


"Το -P (WIN – [0H + PIN = ФЛ (А.5) 


Where N*(u) and № (м) denote the limiting values of the function N(z) when z tend to the 
cut (0,1) from right and left side, respectively. We can rewrite Eq. (A.5) in the form 


G(u)N*(u) — N-(p) = εί) (A.6) 
where 
Q(u) = Ῥίω) 
С) = оиу 7р р) | το 
and 


Df 
eil = OG) + DH) e 


The function N(z) defined by Eq. (A.4) has the following properties: 
a) N(z) is analytic in the complex plane with cut from 0 to 1. 


b) N(z) is vanishing at infinity as ы 
= 


ü а r with &«1l as z>0 and 21, respectively. 
== M 
Now the problem of solving Eq. (A.1) is reduced to the solution of the Hilbert problem 


(A.6) with the additional conditions (a), (b) and (с). 


c) N(z) is bounded by £ and 
2 
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Ц is easily seen that the function X,,(z) defined by the equation 


1 
1 
— I1nG | 
Χμία) = en poe du (A.9) ) 
ch AES 
0 
is a certain solution of the homogeneous Hilbert problem that is of Eq. (A.6) with g(u) =0.. 
Let us investigate the behaviour of the function near the ends of the interval (0.1). . 


It can be easily proved that near the points z=c, where cj —0 and cı = 1, respectively : 


Se: e (= — c, ec? Q(z) (A.10) | 


where 


In Ge, 
жи + 1Ви = (— 1)"+1 SE 


гл: 


(А.11) 


апа Q(z) is bounded and nonvanishing at the point c,. From the definition С(и) for u € (0,1) 


Tee lees № ot E 2i9(u) (A.12) 
where 
Olu) = arg [Ο(μ) — PW] (4.13) 
It gives 
In G(0) — In G(1) — 0 (A.14) 
Hence we obtain that 
αρ = αι = В = В, = 0 (A.15) 


and the index of the singular integral equation (A.1) x is equal to zero 


x#=%+a, = 0 (A.16) 


So we can conclude that А (=) is a bounded function near the points c,. 
Moreover, from the definition of Xp(z) we see that | 


--- 


lim Ae! = constant > 0 | (A.17) 
& +00 Ν 

As the function Х,, (2) is the solution of the homogeneous Hilbert problem it satisfies 
the equation 


Xk 
АШ и) (А.18) 
Ав (и) 
Inserting it into Eq. (A.6) we obtain 
ХЕ — (ЖЖ (ш) = eu) Xx) (A.19) 
Defining a new function 
M(z) = N(z)X,,(z) ~ AI 


e can write Eq. (A.19) in the form 


M*(u) — M-(u) = είω) Ху (и) (A.21) 
he general solution of this equation has a form 


1 


M(z) = „ле y+ P,(z) (A.22) 


y — zZz 


here //2,(z) in an arbitrary polynomial of the n-th degree. Hence we have 


1 
κ ε(νΧπ)., Pala) 
2πι Katz) А y ----ᾱ X1(z) 


(A.23) 


ow we must test whether the function N(z) fulfils the conditions (a), (b) and (с). It can 
e seen immediately that the conditions (a) and (c) are always satisfied. In order to satisfy 
he condition (b) we must put 


(p,(z) Ξ0 (A.24) 
Using Eqs. (A.4), (A.23), (A.24) and Plemely formulas we obtain the solution of the 
ingular integral equation (А.1) in the form 
Qiu) 
Vu) — Du) 
1 
LOW = DWI TEE Ge 
πι Xin (и) J 100) + Фо)? — u] 


Yu) = Kr{f(w} = flu) + 


Now И we put f(u) in Eq. (A.25) in the form 


20 
ϱ 0 a ης 
—с ТУТ C] — Ck ve і Y(v) 
flu) = FE оде du, эйе Deh { “dv 


(A.26) 
we obtain the Fredholm type equation for the function Y(u). 


APPENDIX B 


The method of evaluation of the Fredholm integral equation (2.24) and of the critical 
condition (2,21) 
Let us write the Fredholm equation (2.24) in the form 


1 _ 2p 
Уд — fe KO, μ) γώ» =e CA) Br Ay) (B.1) 


0 
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1 | 
; ' | 
| ота μι S iW (аһ en (В.2)) 
КИДЕ 2T(u) | 9 I — u)(n-+ v) v+u 2) 
0 
1 
А C1 — бле PR bn SA [2.57 — n) 1 B3 | 
A) = pron S(u)d (u, — v) De, Τμ) «EROS 7 ( | ) 
and 
Я 
c [οι — с]? / η (η) Φο(η. νι) fi B. 
Pip) = с, τ᾽ 5009, шь у 2c,T(u) , n—u der 


Eq. (В.1) can be written in the spherical representation. U sing Py-approximation we have: 


ene ааа (В.5) | 
m=0 
where 
Им) = Σ Vm D,(2u — 1) (B.6) | 
fil и) = CG > e D, (2 zu E 1) (B.7) 
N 
(м) = 2) Du (2и — 1) (В.8) 
т= 0 
апа 
12 | 
SW fe 20 
К (0) = (2n+ 1) | | ο WK Kb, uYD,, (2v — NP, (и — Паша 
ὃ ὁ 
1 
20 - 
—=-(2п KD | ο № (ον — 1)а» | K(v, н) (и — Пан We 
0 ô | 
Eq. (B.5) сап be written in the form 
E ie ο 
IA. Ge А 00)! II = | e ub 22+ е Gë e B.10) 
The solution of this equation has a form 
Ke жыр 
μα... ЦО 5 ee € = el SC e В.Ш 
| 
Hence ultimately the function Y (м) 
ος e E 
Yu) = |е get oP er “ο Oa D (12) 


| 
1 
Where μα 


ΙΤ . 
mall denotes the transposed matrix. 
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Finally the critical condition (2.4) has the following form 


l 
| * 9 
E. o e m S 
SERES dv ой А e = 7) f y ‚ pers 
e Lm 3 ganio М9 d Κμπίο)ι - Dm(2v ^ 1) E) dv 
i РИ 
0 
RR zB. 
e lvi p hn 
(В.13) 


The transcedental equation (B.13) involves only material parameters of the system 
onsidered and the thickness о of the multiplying slab. A similar method can be used to 


he solution of Eqs. (3.45), (4.36) and of the critical conditional (3.42) (4.31), respectively. 
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In a previous paper (Acta phys. Polon., 19, 759 (1960)) the author presented purely geome- 
trical hypotheses about the anisotropic propagation of spin waves in regular antiferromagnetic 
superstructures. In the present paper theoretical proofs for the rightness of these hypotheses 
are given. 

The general dispersion law for long spin waves has been derived in a tensor form, starting 
from the Ziman's (Proc. phys. Soc., 65, 540 (1952)) relation based on the Holstein-Primakoff 
approximation. The formulae for all 7 regular superstructures have been obtained, valid for 
arbitrary values J, and J, of exchange integrals between nearest and next nearest neighbours 
respectively. For the extreme cases (J, = 0 or J, = 0), the dispersion anisotropy caused by the 
superstructure appears to be very pronounced, in full agreement with above mentioned hy- 
potheses. For the more realistic non-zero values of exchange integrals, the anisotropy is also 
evident but depends on the ratio J,/J,. More complete discussion of the shape of isoenergetic 
surfaces will be published in the second part of this paper. 

The theoretical treatment leads to the prediction of a new phenomenon in antiferromagnet- 
ies, namely the anisotropic spin wave dispersion which must be attributed to the monoplanar as 
well as biplanar type of regular superstructure. It has been shown that the magnetic anisotropy 
does not affect in principle the phenomenon under consideration. 

Some remarks have been given concerning the possibility of occuring of analogous phe- 


nomena in other lattices and superstructures. 


I. Introduction 


The purpose of the present paper is to verify theoretically previous hypotheses of 
he author concerning possible anisotropic behaviour of spin waves in individual cubic 
egular superstructures. These 10 hypotheses (see also Sec. 2. C. of the present Chapter) 
ауе been formulated on the basis of geometrical considerations in the paper by Cofta 
1960c), which will be called here I. 
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According to the principal thesis of I, т the regular! superstructures (excepting 
natural ones) the energy of the spin wave depends on its direction. The superstructu 
-dependence of spin wave propagation has been foreseen for the first time by the аш 
(Cofta 1957) by reason of approximate calculations carried out for a few superstructu 
only, and disregarding the effects of neighbours more distant than the nearest ones. A note 
ee the influence of the second (i.e. next nearest) neighbours has been publishe« 
some months ago? (Cofta 1960a). | 

In Sec. Г. I the scope of our theory is determined; in Sec. 7. 2 the possible origins o» 
anisotropic propagation of spin waves are shortly written down. The dispersion relatiori 
obtained previously in the harmonic oscillator approximation 13 adapted to our purposes 
in Ch. IT; here the use of power series expansion of the square of energy is justified (Sec. IT. 27 
and general formulas for the anisotropic dispersion tensor are derived (Sec. II. 3). Som 
general conclusions deduced from these formulas are given in Sec. ПГ. 1; the results of 
detailed calculations are presented in Sec. ЛТ. 2 for the case of isotropic superstructures 
and in Sec. 11. 3 for the interesting case of the anisotropic ones. In Sec. JV. 1 all the results: 
are summarized. Some additional remarks are given in Sec. IV. 2. In the Appendix I the¢ 
use of an real function B(k) is justified, in App. П the identity of semiclassical and Ziman’s: 
results is shown, and in App. Ш the applicability of the power series expansion of the energy: 


is discussed. 


[. 1. Restrictions and assumptions 


The subject of our interest are the single domain monocrystals, represented by the: 
ideal antiferromagnetic translation lattices (i.e. Bravais lattices) of equal tons. We consider: 
here only those lattices which are composed of two identical interpenetrating ferromagnetie > 
sublattices, the magnetization of each being antiparallel to the other one. The magnetic mo- - 
ments are assumed to be in practice the spin moments, localized in lattice points and coupled | 
by the isotropic exchange interactions of the form 


A 


#- зуу, 75,5, (0. 
to] 


where hS, is the operator of spin angular momentum at the i-th lattice site, the summation 
being extended over all the pairs of spins. In other words, we suppose the Heisenberg- 
Dirac Hamiltonian, involving the direct as well as all types of indirect exchange interactions 
which may be written in the form (1) (e.g. the superexchange, see Anderson 1950). Let us ᾿ 
call all so generalized coefficients Jj "exchange integrals" and restrict our considerations ` 
only to the isotropic ones. This means that we suppose an appropriate high symmetry of 
the chemical lattice as well as that of the ionic orbitals, excluding all possible causes of 
directional dependence of the exchange integrals. Moreover, we shall consider solely the 


* - 


1 As to the definition, see Cofta (1960c). 
2 Moreover, a short review concerning th 


e problem of superstructural anisotropy of spin-wave propagation 
has been published recently (Cofta 1960b). 


D 
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nearest and the second interactions (i.e. the interactions between the nearest or next nearest 
E cighbours respectively), although the formalism used here permits to include the effects 
ES arbitrarily distant neighbours. We shall neglect the influence of crystal boundaries. 
4 We shall also assume the g-factor to be isotropic and equal for both sublattices. Our 
“considerations will be confined to the cubic lattices only (see See. Г. 2. А) with regular 
A in superstructure (see Sec. /. 2. C). Our treatment will be valid in principle for tempera- 
tures which are sufficiently low as compared with the Néel point (see Ch. 77) 


I. 2. Possible origins of anisotropic propagation of spin waves 


NOW VUA és tz vw 


Obviously, the superstructure is no unique possible origin of anisotropic propagation 


of the spin waves; other causes may also exist, compensating the considered superstructural 


Ü 
effects. Let us comment shortly upon all these origins, as listed below. 


wed 


А. Anisotropic structure of magnetic lattice 


The anisotropic structure of magnetic lattice presents a quite evident cause of aniso- 


ΑΝ ΑΛ 9 


tropic behaviour of all kinds of the wave motion through the crystal. In the case of spin 


4 


"waves, the anisotropic structure causes the directional effects even in ferromagnetics. In 


order to eliminate the structural anisotropy, we shall confine ourselves to cubic magnetic 


structures only. 


B. Anisotropic interactions between atomic spins 


АТА 


Many mechanisms may occur, causing this sort of anisotropy. Let us survey below the ` 


essential cases of such a type. 


a) Multipolar and pseudomultipolar coupling 


NY Υπ 


'The dipole-dipole as well as quadrupole-quadrupole interactions caused by the spin- 
orbit coupling may modify slightly the dispersion law for spin waves even in the ferro- 
magnetic case (Kowalewski 1960). But the anisotropy of that origin amounts to only about 
2 percent of the isotropic part of the dispersion law. Supposing that the analogous effects 
modify the antiferromagnetic dispersion law not more than in the ferromagnetic case, we 


- shall ignore the multipolar interactions in our considerations and assume the Hamiltonian 


of the form (1). 


b) Anisotropic exchange integrals 
N 


The exchange integrals may depend on the direction on account of several different 


causes, namely: 
a) when the chemical lattice has lower symmetry than the magnetic one; such a situation 
may give rise to different indirect interactions at the same distance (a simple two-dimensional 


example is shown at the Fig. 1) 
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В) when asymmetric orbitals of magnetic ions are present; the coplanar orbitals (Goode- : 


nough and Loeb 1955) may be here a typical example; 
y) when different sublattices are formed by different ions. 


οί )eC)o 
e( )o()e 


o()e()o 


Fig. 1. Example of a chemical lattice in which anisotropy of exchange integral can occur. Large circles symbolize 
the anions, small — the cations with opposite spins (white and black respectively) 


All the mechanisms listed above may cause anisotropy of exchange integrals theoreti- 
cally even in the ferromagnetic case. The assumptions made in the present paper, eliminate 
all causes of this art. They will be the subject of future investigation. 


C. Anisotropic superstructure 


As we have pointed in the paper I, all regular superstructures excepting the natural 
ones have an anisotropic character. Our restrictions and assumptions admit solely this cause 
of anisotropy so that all anisotropic effects obtained in our calculations, are only of purely 
superstructural origin. Let us recall shortly our previous presumptions concerning these 
effects, presenting them in a tabulated form. 


The expected character of the 
Lattice and superstructure - - —— - === 
nearest interactions term second interactions term 

sc— natural isotropic isotropic 

ze (001) Е (P) anisotropic [001] (P,) anisotropie [001] 
se—(110)—(110) (Рз) anisotropic [001] (P4) anisotropie [001] 
bce— natural isotropic isotropic 

bcc— (110) (P,) anisotropic (P) anisotropic [001] 

[110] and two other directions 
сс — (001) (P,) anisotropic [001] (P,) isotropic 
fcc — (111) (Py) anisotropic [111] (Pio) isotropic 
Explanation: 


We use the notation of regular superstructures introduced-by the author in the paper [. To characterize 
the anisotropy we have given in each case the preferred direction which should be expected (mentioning 
of one direction only means the uniaxial anisotropy). The symbols P, number the hypotheses in our paper I. 
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The above hypotheses concern only extreme cases i.e. the case of vanishing second 
neighbours exchange integral (nearest interactions term) or that of vanishing nearest neigh- 
bours exchange integral (second interactions term). Suppositions as to the form of the disper- 
sion law for J, #0 and J„#0 simultaneously have not been possible, considering that the 
resulting effect is no simple superposition of the two special cases mentioned above. 


П. Adaptation of the dispersion law 
II. 1. The approach used 


| The interactions between spin waves play practically no role in low temperatures, 
probably up to the 3T. We will develop our theory for this range of temperatures (see also 
Sec. IV. 2). In attempt to investigate only the propagation phenomena of spin waves without 
their interactions one another, we have to employ, of course, the harmonic oscillator approxi- 
mation. There exist several different methods using this approximation: 1° the approach 
of Holstein and Primakoff, 2? the approach of Dyson (both in the so called “linear” or “free 
magnons” approximation) and 3° the semiclassical approach. All three methods mentioned 
above lead to the results which are exactly identical to one another, which was to be expected 
considering the identity of approximation?. Therefore, each of these three approaches can be 
equally chosen for our purpose. There are available formulae obtained by the method 1° (Zi- 
man 1952) or by the method 3° (Cofta 1959). Both formulae differ merely in notation. Let 


us employ here the elegant notation by Ziman. 


II. 2. Expansion of Ziman’s formula 


We shall investigate the phenomenon under consideration only in the absence of applied 
field. Then, the dispersion law for spin waves in antiferromagnetic regular superstructures 


has the following form in the notation of Ziman (1952): 


E(k) = V(4())? — |BR)P (2) 
where 
A(k) = gBH, +25 > Jfo,„—-%(1+o,) cos kr} 
* B(k) = — S Σ) Διί] — σι) exp (ikr) (3) 


3 The above statement holds for ferro-, antiferro- and ferrimagnetics. The identity of results obtained 


by the methods 1° and 3° for the case being treated here, is shortly discussed in the Appendix II to the present 


paper. The proof for the more general cases will be published in the subsequent papers. It is to be noted that 


(for the case of antiferromagnetic natural superstructures with nearest 
Kubo in the paper by Nagamiya et al. (1955). The identity of results 
ins (antiferrimagnetism) has been obtained 


the first remark concerning this identity 
interactions only) has been published by 
of 1? and 2? for the ferrimagnetic simple lattice with antiparallel sp 
by J. Morkowski (private communication). This is in agreement with the result of J. Kociüski (1960) for the 


ferrimagnetic case with parallel spins. 
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Here 5 denotes the maximal value of the z-component of atomic spin, J, — the exchange 
integral between the arbitrary spin and its neighbour in the relative position т. Further 
more — the / is the energy and k the propagation vector of the spin wave; f stands fo» 
the Bohr's magneton and g for the Landé's factor. The number c, has the following property i 
0,=1 if the spin at the point r belongs to the same sublattice as the given spin at the poini 
r—0 ; otherwise g,— —1. (That is the Luttinger-Ziman method of description of the 
regular superstructures). The effects of magnetic anisotropy are taken here into account 
in the phenomenological manner, H, being the absolute value of the effective magnetic 
anisotropy field, directed along the z-axis. 

The Ziman’s formula holds for all regular superstructures (Cofta 1960c). By reason 
of our assumptions, the function B(K) is a real one (see Appendix I) and we can drop the 
symbols of absolute value in (2). | 

In the temperatures for which the free magnon approximation is valid, only the spin 
waves with small values of k play the essential role: 


(a being the lattice constant). Therefore, we apply for our purpose an approximate form! 
of the formula (2). Using the Taylor's expansion? 


dF SW el 
k 


5 (k © k) $ 


FOS FQ) κ. | 6) 


dk | „о dk?) 0 
for the function E, we obtain 
E? = Ck- Mk 4... (6) 
where the constant C denotes 
C = еВН (ep Ημ + 2B) 
B, = Β(0) = — S$ J, — в.) (7) 
T 
whereas the tensor М may be expressed in the form 
M=SNJFlror){Al+o)-B(l- с,)! d (8) | 
= 
with 
Ay = gBH, + В, г (8) 
It can be easily determined that the third non-vanishing term in (6) i.e. the А term amounts ! 
аё? | 


to ds times the second one. Since on account of (4) 


a?k? 


15 ΣΠ 


4 То avoid mistakes, we employ here the symbol © for the dyadic multiplication of two vectors, and the 


dot for the scalar multiplication of two tensors of any rank (order). Note, that for anv second rank tensor T' and 
any vector e the relation (v © в): T = п Tv holds. á | 
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we can ignore all higher terms in (6) and write in a very good approximation 


E ~ {C+ ΚΜ) 


wie 


(9) 
3 Sometimes the dispersion law for long spin waves in antiferromagnetics is expressed 
also in another form, namely that arising of the approximate power series expansion of the 
mergy. In our general case we get thus: 


АМЕ 


E(k) ~ VC A = 
(k) ~= y C (10) 


The shape of isoenergetic surfaces is, of course, the same for both the forms. Nevertheless, 
e prefer the form (9) because it is valid without any restrictions as to the values of Ё except 
of (4). Due to this fact, this form is convenient for simplified calculations, neglecting the 
jagnetic anisotropy. The form (10), instead, behaves differently, since the infinite Taylor 
series of E has a limited convergence radius?, depending on the magnetic anisotropy field 
strenght (see Appendix III) and, moreover, the approximation (10) (ie. the dropping of 
"the higher terms) is applicable in a rather more restricted range namely for ak<10=? pro- 
vided the usual values of Ἡ (see also the Appendix III); for example, in very thin antiferro- 
‘magnetic films this condition is not satisfied. For vanishing magnetic anisotropy (C=0) , 


E form (10) is obviously no longer valid. 


II. 3. The anisotropy tensor M 
We see from eq. (9) that the anisotropy of the spin wave dispersion is characterized 
completely by the tensor M. We will now adapt the formula (8) for the case to be investi- 
‘gated here. Let us denote by J, the exchange integral between the nearest neighbours, 


“by Л however that between the second ones. Then, by virtue of assumptions of the Sec- 
tion I. 1, we may write: 

By = — 250,2, Ла, ) (11) 
E z, and z, denote the numbers of nearest and second neighbours respectively, having 
the spins aligned oppositely to the spin at the point r=0. Moreover, let us introduce the 


_ notation 


£4 = gpl, (12) 
From (11) and (8) we get 
i M Уро) CT Or) (13) 


where we have used the abbreviation 
| G, = ад ο. 25By) о, (14) 
о 


5 or, strictly speaking, the limited convergence ellipsoid. 
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n 5 

The symbol >J in (13) is used for the summation over the nearest neighbours only, and ^l 
over the second ones. 

In order to verify our hypotheses presented in the paper I et us divise the tensor M(J, , J) 


which is a function of Л, and Л, in three parts: 
М = M" + M+M” (15) 
defined as follows: 
М" = M(J,, 0) 
M' MOD (16) 
M"—-M—M"—M* 
Now, by means of (16) we may treat the k-dependent part of E? as separated into three 


terms: the “nearest interactions term” kM”k, the “second interactions term" АМК and 


the “mixed” опе: kM”k. From (13), (14) and (16) we have: 


M" = J,‘ >) (ror) {e,(1+0,) — 45°/ z7o,) (17a) 
M: = J, ron {= (1+0,) — 453/270) (17b) 
М" = — 4527 7,2 >) (r Orjo,+ 27 2> roro} (17e) 


Ш. Results 
III. 1. General conclusions . 


Two general conclusions may be deduced from our formulas. The first one concerns 
the effects of magnetic anisotropy on the phenomenon investigated here. Coming back to 
the formula (14), we see that the first part is always isotropic (see also the Section Л. zu 
In the second one, however, the energy of magnetic anisotropy oceurs only as an item of 
the sum ед +2SB,. The magnetic anisotropy energy is known to be roughly 100 times smaller® 
than the exchange energy ep (defined below): 


вд 77 10736, (18) 


For the case of non-vanishing second. interactions we have to define the “exchange energy" 


in a rather more general sense as usually, sav: 


ёв = 45? { |, e τ dëse (18) 


For most antiferromagnetic substances holds: 25В уер so that, generally speaking, in the 
second part of (14) we can neglect the £g in comparison with the exchange term 2SB,. Hence, 
we can conclude that the magnetic anisotropy can play an important role in the investigated 


* E. в. for the MnF, we get £472 X 10? x 258, ~2 х 1072 eg, see Johnson and Nethercot (1959). 
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superstructural effects only in the rather exceptional case 25B, «e4 which implies on account 


of (11), (18) and (18): 
Juss + Jar S es + Wiles} 107 19 


In that case e, can become the unique factor in both the parts of the formula (14), namely 
when 


B, w0 (20) 


We will show in the Part II of this paper for each superstructure separately that (20) is 
rather not to be expected. 

The second conclusion concerns the behaviour of the E(k) when the coordinate system 
is changed. Let us consider a given lattice with a given superstructure, which is oriented 
in an other direction than that quoted in the table 1. Further let us apply the formula (9) 
to the rotated lattice. It can be seen that the only difference in the dispersion law in compa- 
rison with that for the unrotated lattice is the change of the orientation of quadric surface 
of tensor M, the new orientation being consistent with that of the lattice; the essential 
parameters, however, as well as the constant C remain unchanged. In other words, the disper- 
sion relation does not depend on the relative orientation of spins (and the field H4) with 
respect to the whole lattice. Hence, the SA, as well as 5? involved in (9) may be treated 
as the invariants of rotation of coordinate axes. This is obvious when we consider that the 
approximations used here imply small deviations of spins from their common direction, 
so that SH, and S? behave as the scalar products of two vectors. 


III. 2. The isotropic dispersion 


Let us calculate at first the tensor M for isotropic superstructures i.e. the natural ones. 
We have, of course, in this case z7 =z, (2, being the number of all nearest neighbours) 


and z; =0. 
Then from (14) and (11) we find: 


G = ro, E 52,0, (21) 
The definition of natural superstructures implies also o, — —1 for all nearest neighbours 


and o,=1 for all the second ones; then the formula (13) yields: 


M = 45272, У) (r o v) + J (284 — 45*J,2,) той (22) 
It is easy to evaluate the sum Zur ОТ) for neighbours lying on a sphere with radius r: 
Dur = izr*E. (23) 


Here z is the number of neighbours lying on the sphere and E denotes the unit tensor. 


Therefore, for the case being considered we find 


k -Mk = Mk? EEN 
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with 


М = 2527222 t+ X Ja (264 — 4S2], )r? (25) 


n?n'n | 3 n^n 
where r, and г, denote the radii of nearest and second neighbours spheres respectively, Ze 
being the number of second neighbours. The result (25) shows tbat the magnetic anisotropy 
cannot play decisive role in isotropic dispersion relation, since е. is always much smaller 
than S?|J,|z,. In practice, we can disregard е, in expressions for М as well as in that for 


the constant C: 


C = gH (ЕВНА + 280) ~ 22BH,B, (26) 
because for natural superstructures 
В: = — 2572, (27) 


and therefore By never vanishes. Hence, for the case of natural superstructures of cubic 


lattices, the spin waves obey the isotropic dispersion law: 


E=VC+ Mi? (28) 
where 
C= — de uz, : 
Mise Is) (29) 
If |J, |J, |. then the relation (28) can be written approximately in the form: 
E22 | rs V zin 3 ы rok (30) 


TL Anisotropic results 


We will present now the results of calculations of 6 components of the symmetric 
tensor M individually for each of 5 anis« tropic superstructures of cubic lattices. The lattice 
constant will be denoted by a. The figures annexed to the present section show the relative 
positions of the nearest and second neighbours’. 

a) The case of sc-(001) 

From the Fig. 2 we find: z, =2, ᾱ- =8; hence Ву = — 45(7, +47). From the 
formula (17a) we get: 


М» = My, = AaY(e, — 48°J,)J, 


Μὴ = 1602527? (31) 
Mey = My, = М», =0 


Considering (18) we find a pronounced uniaxial anisotropy of the nearest interactions term, 
which may be written in a good approximation as proportional to 


ky t kè — 13 (32) 


+ — 


? The more complete illustrations of superstructures may be found in the paper Г. 


E 
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Since the preferred direction is the [001] one, our hypothesis P, is hereby fully vindicated 
theoretically. | | 


Likewise we obtain for the second interactions term: 
$; о 
Μ᾽, = MES = Ba?e,), 
MS, = (16a) 22J (33) 


М = М, Mag Ὁ 


Since MS, << М°., the direction [001] is very strongly preferred in this case (see also the case 
of fcc-(001) in the section 2d of the present Chapter). In this manner the rightness of our 
hypothesis P, is also completely demonstrated. | 


2 


о 


Fig. 2. Relative positions of neighbours in the case of superstructure sc-(001) 


Taking into account also the mixed term, we get the complete result: 
fus : 22 
M,, = М, = 4a 9, τ. АВ 47,)} + (34) 
ος 2 
И. 16225*(J op AJ) 
the non-diagonal components being equal to zero. This result shows also the uniaxial aniso- 
‚ору with the preferred [001] direction. 
b) The case of sc-(110)—(110) 
According to Fig. 3 we have: z, =4 2, = 8, hence By = — 81 SST For the 
diagonal components of the nearest interactions term we get 


Mr, = М", = 322252]? 
М = да? {ел — SSI, 


(35) 
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The non-diagonal components vanish. The uniaxial anisotropy. privileging the direction: 
[001] is here evident, in full accord with our prediction P}. 


The second interactions tensor is also diagonal in this case, with the same direction 


preferred : 


M? = М}, = 8a%e,J, 
М}, = 256a?82]? 


(36) ) 


Similarly with the result (33), here also a strongly pronounced anisotropy occurs, which ı 


argues for the reality of our hypothesis P,. 


Fig. 3. Relative positions of neighbours in the case of superstructure se-(110)—(110) 


The full result shows the same type of anisotropy as both the special cases given above: 


М... = M,, = 8a%{e,J,+ 4527 (J, + 2J,)} 


xx Уу 
Мы = Ча (e4J, + 8S*(4J, — JAJ, + 27,)) (37) 
M, = M,, —M,,-0 : 


Let us emphasize, that for some values of ЛАЛ, the effects of nearest and second interactions 
may be mutually compensated ; e.g. when e, is disregarded, we obtain for J, —2J, the isotropy 
of spin wave propagation: Ma = М, 
c) The case of bee-(110) 
As S : ` <. ы E x E „— CS aj "a 
8 shown ut the Fig. 4, we have here Za =, =4, hence By = — 8S5(J, +J). For 
the case of nearest interactions only we find: 
A Nn v € 
Мих ΜΕ = ΜΗ = 2a*e MJ, 
n^ doa 
My, = 2a*(e, — 1682] )J, (38) 
AR — TE I | 
Му = М" =0 
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| Fig. 4. Relative positions of neighbours in the case of superstructure bcc-(110) 


In order to show more clearly the anisotropy appearing here, we have to go over to the new 
earthesian coordinates axes x’, γ΄, 2’ in directions: [110]. [110], and [001] respectively. The 


appropriate transformation is given by 


К =Tk (39а) 
where 
wel) 
о (39b) 
ү GE 


In the new coordinates the componente of M are the same аз the components of а tensor 
MW = TMT- in the original system, namely 


(MZ)! = Mz, — Mg, = 829259, 

(Mz) = М" MZ, = 4a?) (e4 — 85°J,) 

(Mia) = Mz, = Zeie, 

(Мы) = (Mi) = (М) =0 E 


The result (40) gives proof for the anisotropic behaviour of the nearest interactions term, 
one of the three favoured directions being perpendicular to the ferromagnetic planes, in full 


accord with our prediction P;. 
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For the second interactions term we get in the new coordinate system (39): 
p 7 / DI aD p! 
(Mz) = (My) = 32a?S?J, (41) 
(M;,)' = 4a7J,(e4 — 8527.) 


the other components being equal to zero. The evident uniaxial anisotropy in the [001] 
direction, which may be seen from (41), gives the theoretical verification of our hypothesis Pg. 


. . » . 
When, moreover the mixed term is taken into account, then the full tensor M appears: 


also diagonal in the system (39): 
Р ИГ ЭА 
М,, = 32a2S2(J? — J?) + 4a%e,J, (42) 
M,,— — 32a?S2J (J, + J) + 2a7e4(J,, + 2J,) 
d) The case of fcc-(001) 


From the Fig. 5 we get: 2, —8,2; —0; hence B, = — 1657,. The formula for the: 


nearest interactions term gives following components of diagonal tensor M”: 

Je πα E 

М». = MS = 2a*e,J, 

Ме = Gaeta po (43) | 
As may be seen, this result is in principle identical with both the ones given by (33) and (36). . 


This is no surprise, since the substructure of sc-(001) as well as those of sc-(110)—(110) | 


Fig. 5. Relative positions of neighbours in the case of superstructure fcc-(001) | 
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form simply fcc lattices with the regular spin ordering (001). The result (43) demonstrates 
the rightness of the hypothesis P,. On the other hand the essential role of the magnetic 
anisotropy in the components M,, and M, in (43), (33) and (36), contrasting with its role 
‘in the isotropic case, seems rather unexpected. (The result (40) shows also a similar beha- 


viour of AM A 


For the second interactions term our formulae give isotropic result: 


Mi, = M3, = Mi, = 4a%e,J, (44) 


M3, = Mi, = Ms, =0 


which corroborates the prediction, formulated in our hypothesis Pg. It is worth while to 
note the quite exceptional role of the magnetic anisotropy in (44), which seems to suggest that 
the exchange interactions of the form (1) between the second neighbours can never be the 
unique coupling mechanism in the case of superstructure under consideration. Only the 
simultaneous presence of the nearest and second interactions gives rise to appearing of 
purely exchange terms in components of the tensor M (see the result (45) below). 

The full expression for M has the form 


M,, = M,, = 2a%e4(J,, + 2J) — 64α553],]. 


M,, = 4a?e 4J, + 644352] (J, — I) (45) 
nope В M 0 


Fig. 6. Relative positions of neighbours in the case of superstructure fee-(111) 
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e) The case of fee-(111) 
In this most interesting case of superstructure we have 2, =z, =6 which implies 


By= — 125(7, + J;). The calculations of the sums (17) (see the Fig. 6) give for the nearest 


interactions term: 


Ми, = М”, = ΜΗ = 2а2е 4], ad 
М», = М", = М", = aJ (e4 + 2453] ) 


In order to show the anisotropic nature of the non-diagonal components we introduce another 
coordinate axes, the new axis z having the [111] direction in the original system. The two 
other axes may be chosen quite arbitrarily (provided that they form an orthogonal system), 
because the anisotropy being considered appears to be an uniaxial one. Let us put the axes x 


and y in the direction [112] and [110] respectively. Then, using the transformation 


να ste aae ya 
τε WE E уЗ (47) 
KATETE γ2 
we get in the new system: 
(My) = (My) = MS. My, = a?J (e, — 2452] ) 
Ж е = М", ч 2M; = 4a*] (e, + 1252] ) (48) 


(Mz,)’ = (Mz) = (M*y =0 


This result demonstrates that our predictions, formulated as the hypothesis Pg, are right. 
For the second interactions term we σοι: 


M;, = Mi, = Мз, = 48а252]2 (49) 


yy zz 
other components being equal to zero. This isotropic result is the theoretical argument for 
the reasonabless of our hypothesis P. 


The full formulae for the case under consideration are: 


Myx = Му, = аєа], + 245%, + JAJ, — J} (50) 
M,, = 4a*1e4J, + 1252(7, + T 


The isotropic result obtained by Elliot and Lowde (1955) is due to an error? in their caleula- 
tions, 


IV. Conclusions and remarks 
IV. 1. Summary of results 


1° In the paper I we have presented on basis of geometrical considerations the 10 hy- 
potheses (see the table in Sec. 7. 3. C) concerning possible behaviour of spin waves in the 
1 1Y ‹ ER ra aT ο νο γη] rac far tha r aytr ас 
individual 5 regular superstructures for the two extreme cases (J, = 0 ος J, =0). In the 


This isotropy is a consc quence of the faulty assumption that in the factor A—|B| the sum over nearest 
neighbours vanishes, І 
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section III. 3 of the present paper, we have given the theoretical proofs of the rightness of 
all these hypotheses in all their details. The superstructural anisotropy of spin wave pro- 
pagation predicted in P,, Р,, P4, Ра, Ρε, Pg, P; and Р, appears very pronouncedly. The 
preferred directions agree in all cases with these expected. Also the uniaxial character of 
the anisotropy is shown for all cases when predicted (i.e. with the exception of P,). The 
isotropic dispersion relations expected by Pg and P4, have been proved. 

2? [t has been supposed that the anisotropy of spin wave propagation will be practically 
quenched by the energy of magnetic anisotropy. The theoretical calculations presented in 
the Chapter ПТ remove entirely such apprehensions. The existence of magnetic anisotropy 
does not affect the superstructural phenomena being considered here, except in the case when 

AS? Spa Js, | 5 24 (51) 
which seems to be rather impossible for all regular superstructures (a closer examination 
will be given in the Part II of the present paper). Only in this case the dispersion relation 
might be the isotropic one, independently of the superstructure. 

3° From the formulas obtained in ΠΠ. 3 one can get easily all the partial results published 
in our previous notes, quoted in the Ch. 7. This holds also for the results obtained by means 
of the kinematical method (Cofta 1957)? 

4° From the result (44) we can conclude that in the case of /сс-(001), the second in- 


teractions play probably an unimportant role. 


IV. 2. Final remarks 


It is to be noted that the extreme cases investigated here as nearest an second interactions 
terms respectively, seem to be rather unrealistic. We have calculated these terms solely 
in order to justify our hypotheses of the paper I: 

In practice, neither J, nor J, vanishes, and the resulting expression involves both the 
corresponding terms as well as the mixed one, which has in general a rather complicated 
form. In consequence, a compensation of nearest and second interactions effeets can occur 
for definite values of J, and J,, so that the dispersion relation can become isotropic. Such 
an effect is possible even when both the interactions favour the same direction as we have 
seen in IIL. 2. b. We shall analyse accurately the shape of isoenergetic surfaces in the Part II 
of the present paper, to be published shortly in the A.P.P. 

It is to be suggested that also the anharmonic terms in the Hamiltonian must show 
a kind of the direction-dependence arising from the superstructure. For the same reasons 
as those presented in the paper I, these higher order directional effects must show, of course, 
an analogous character (i.e. favouring the same directions) as the ones investigated here. 
On this ground one can hope that the phenomenon of the superstructural anisotropy may 
be observed in a larger region of temperatures than that one limited by the free magnons 
approximation. Moreover, some directional dependence of the resonance line width can be 


expected for the same reasons. 


9 As to the equivalence of the kinematical method of Keffer, Kaplan, and Yafet with the other semiclas- 


sical ones, see the paper by Cofta (1959). 
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At the end, let us note some possibilities of extending of present investigations. First ! 
of all, the analogous effects i.e. the superstructural dependence of the spin wave propaga- : 
tion, must occur also in the cases of anisotropic interactions. Especially, the case of the: 
anisotropic exchange integrals is an typical domain for such investigations (some results : 
will be given shortly). Further, the superstructural anisotropy of the spin wave dispersion 
is also possible in the non-translational antiferromagnetic lattices (i.e. in lattices with a basis). . 
Here the effects will be probably weaker than in translational lattices. 

Effects similar to these obtained in the present paper will appear also in some kind on 
non-symmetric superstructures, namely in those in which the ferromagnetic planes exist 
(excepting the isotropic cases, analogous to the natural cases of the regular superstructure). 
The approximate calculations for the nearest interactions only lead to the strongest aniso- . 
tropy possible: the direction perpendicular to the ferromagnetic planes is a “forbidden” 
one (see also Cofta 1959). Further investigations seem to be very interesting. 

The superstructural phenomena must occur obviously also in ferrimagnetic lattices. 
On account of overlap effects, the character of propagation anisotropy is here surely different 
from that in antiferromagnetics. The case of ferrimagnetic non-translational lattices is, how- 
ever, interesting in historical respects. In this most complicated case, the anisotropy of 
dispersion relation caused by the superstructure has been obtained for the first time, namely 
by Kouvel (1955) in his excellent spin-wave analysis of magnetite lattice!?, Here the super- 
structure in B-positions causes a weak anisotropy in the [110] direction. If the usual estima- 
tions of exchange integrals in spinels are made, this anisotropy amounts to less than 2 percent 
and therefore seems to be unobservable. 
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APPENDIX I 
The reality of the function B(k) 


By virtue of the assumed isotropy of the exchange integrals, we have 


L,-J, (АТ) 


‘© This paper was not known to the author before 1958 
published as ye 


; this classie work has been, unfortunately, not 
t. The kindness of Dr Kouvel in sending to the author a copy of his work is greatly appreciated. 
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Moreover, since the regular superstructures are symmetric, it follows immediately from 
their definition: 


/ Op=Op (A212) 
Then from (4) we get: 
ES TSS ЛА — орде" УЛ (1-0) *} = —S У, Л —0,) cos kr. (A L3) 


which is а real function. 


APPENDIX II 
The identity of semiclassical and Ziman’s results 
The semiclassical treatment of spin waves in regular antiferromagnetic superstructures 
gives for isotropic exchange integrals (Cofta 1959): 
E — Wan, —») + ph)? — 51/9} (A ILI) 
where 


fk) = 2:3 J, (cos ler — 1) 


T 


f(k) = 2 > J,cos kr (A 11.2) 


о. 


+ E 
and y — 0 denotes the magnetomechanical ratio. The symbols >) and >) are used for the 
summation over the neighbours respectively parallel or antiparallel to the given spin (Such 
a notation is very useful for practical evaluation of sums over the neighbours). In the Ziman’s 
notation, the sums (А П.2) may be written as: 


Ak) = Ὁ (1+ oJ, (cos kr — 1) 


T 


АК) = Σ (1 — oJ, cos kr (A 11.3) 


T 


^ => (о =, 


Considering that yh = — gf and taking into account ‘the formula (A 1.3), we get: 
SA) — x] + yh, = — Alk) (AILA) 
Sf(k) = — B(k) 


In this way the identity under consideration is demonstrated. 
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APPENDIX Ш 


The power series expansion of the energy 


In order to examine the validity of the expansion (10) let us consider first the natural! 
‚ superstructures of cubic lattices. According to the eq. (24) we are dealing in that case with) 


a one variable only and the Taylor series has the form 


E= yc M SES (А Ш.) 
m=0 
where 
CH 
bo = 1 
d „(п — 3)!! 
κ πο τρις n — 1, 2, 3, ... 
bo, = 0 


Let us determine the convergence radius Е. of (A IIL.1). By means of the well-known 


= и i (A. Ш.2) 


(A. 11.3) 


methods we get 


From (29) we get 


| 
1 | S? REOR e s ES 


Only the order of magnitude of the k, can be estimated, of course. We have to take into 
account that in the case of natural spin ordering, the exchange integral J, must be inevi- 
tably the negative one. As to the integral Л, a twofold situation can occur. Either LI < || 
and then the sign of the J, is unimportant, or || r~ |J,, and then a simple examination of 
the semiclassical Hamiltonian leads to J, > 0. In both the cases the denominator in (A 111.3) 
is of the same order as |7 |2,12. Since for most antiferromagnetie substances holds (18), 
we can write: : 


e (A 111.4) 


n 


[οἱ 3x10? 
r 
Hence, we gel 


k, ~ (6r, )-1 (A 11.5) 


We have to examine also the question of validity of the approximation (10) i.e. of the 
neglecting of higher terms in the series (A Ш.1). The ratio 


ND 
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E 4 


The third term may be ignored only if it amounts at most to a few percent of the second 
one. That implies the following requirement for the applicability of the approximation (10): 
IMEG 


а» οι Е 
вт |1 7 το (A 111.6) 


This inequality restricts the use of (10) rather more than the condition of convergence. 
In the case of anisotropic superstructures the situation is a similar one, as can be seen 
from the formulas (34), (37), (42), (45) and (50), for a fixed direction. E.g. for sc-(001), 


making the assumption 


Ес гову x 10-* 
analogous to (18), we obtain from (34): 


G 4 x 10-2 


M, z e 


and (provided that |J,|— |J]: 
| DAN συ” 
SEL LU E. 


Max а 


Only if one of the exchange integrals J, or J, is equal to zero (extreme cases), the ratio |C/ M| 
may be sometimes 10? times larger than in (A Ш.4), at least for certain directions. So for 


transversal directions we find from (33) as well as from (36): 


IM = фа? (| 


for the second interactions term and from (43) we get the same equality for the nearest 
interactions term. From (40) we find М” = sol | for the longitudinal direction. In (44) we 
get the same relation for the second interactions term for all directions. Nevertheless, all 
these examples do not remove our objections concerning the correctness of (10), because 
the cases of J,=0 or J,=O in anisotropic superstructures seem impossible in practice. 
Consequently, the form (10) of the dispersion relation may be employed only for k being 
sufficiently small. The allowed values of k depend on ε 43.58. As follows from (A Ш.5) and 
(A IIL6), k must fulfil the condition 


ak < 107? 


if the usual values 25Ву— eg are taken. 


REFERENCES 


Anderson, P. W., Phys. Rev., 79, 350 (1950). 
Cofta, H., Acta phys. Polon., 16, 481 (1957). 
Cofta, H., Acta phys. Polon., 18, 215 (1959). 
Cofta, H., Acta phys. Polon., 19, 405 (1960a). 


612 


Cofta, H., Bull. Acad. Polon. Sci. Cl. Ш, 8, 547 (1960b). 

Cofta, H., Acta phys. Polon., 19, 759 (1960c). 

Elliot, R. J. and Lowde, В. D., Proc. Roy. Soc., A 930, 46 (1955). 

Goodenough, J. B. and Loeb, A. L., Phys. Rev., 98, 564 (1955). 

Johnson, F. M. and Nethercot, A. H., Phys. Rev, 114, 705 (1959). 

Kocinski, J., Acta phys. Polon., 19, 565 (1960). 

Kouvel, J. S., Technical Report No 210, Harvard University, Cruft Laboratory, Feb. 1955. 
Kowalewski, L., Acta phys. Polon., 19, 59 (1960). 

Nagamiya, T., Yosida, K. and Kubo, R., Advances in Physics (Suppl. of Phil. Mag.), 4, 1 (1955). 
Ziman, J. M., Proc. Phys. Soc., 65, 540 (1952). 


Го]. ХХ (1961) А ACTA PHYSICA POLONICA ΤΟΙ 


LETTERS TO THE EDITOR 


HUCKEL’S THEORY: DEPENDENCE OF THE COULOMB INTEGRALS 
ON THE ATOMIC CHARGES 
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Formula (1) is assumed. е is determined so as to yield the correct dipole moment for 
fulvene. 

When applying Hückel's theory to alternant hydrocarbons, all core Coulomb integrals 
(a) are assumed equal. This assumption seems justified because all the carbon atoms bear 
equal charges (Coulson Longuet-Higgins 1947). For nonalternant hydrocarbons, however, 
the charge distribution is more complicated. Different carbon atoms bear different charges 
and, accordingly, the Coulomb integrals should be assumed different. No theory giving 
the dependence of the Coulomb integrals on the charges exists as yet. Wheland and Mann 
(1949), and Berthier and Pullman (1948, 1949), proposed x —4 , where 4 is the z-electronic 
charge on the carbon atom; their choice, however, was quite arbitrary (Pullman Pullman 
1952). The only point that is clear is that а should be increased whenever the electronegativity 
of its carbon atom increases. If this is not done, the nonuniformity of the electronic distri- 
bution is overestimated. Thus, for instance, the conventional L.C.A.O. theory yields for 
the dipole momen of fulvene and azulene 4.8D and 6.4D (Debye units), respectively, whereas 
the experimental results are. 1.2D and 1.0D (Pullman, Pullman 1952). Similar effects are 
to be expected for alternant hydrocarbons when ionized, as then electronic distribution 
is no longer uniform. 

In this note the L.C.A.O. theory with selfconsistent values of both Coulomb (а) and 
resonance (В) integrals 15 used for computing the dipole moment of fulvene. The dependence 
of the core Coulomb integrals on the atomic charges is assumed similar to that proposed 


by Wheland and Mann: 
& = Eq (1) 


where в is a constant. It is our aim to find the value of 8 yielding the correct dipole moment. 
The problem of whether formula (1) generalizes correctly Hiickel’s theory to nonalternant 


(613) 
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hydrocarbons remains open, but were this the case indeed, the value of & as found for fulven: 
should be generally applicable. The resonance integrals are assumed to depend on th: 


lengths and orders of the bonds, according to the formula: 
В = Boll — 2.683(R — 1.397 А) + 3.599(R — 1.397 А)?] (2: 
Е — 1.397 А = 0.1802 — p) (3) 


where p is the bond order, and R the bond length. The constants are those found by Coulsor: 
and Golebiewski (in print) using Longuet-Higgins and Salem's (1959) method. Usually 
an exponential formula us used instead of eq. (2), but the derivation (from the vibrational 
spectrum of benzene) gives no information on the terms beyond the quadratic one. Numeri; 
cally the two formulae are equivalent if the accuracy does not exceed 0.001 А. Formula (2) 
is here preferred as yielding the force constant of ethylene more satisfactorily. (For the 
method of computation, see the paper by Longuet-Higgins and Salem). 

The numerical results are given in Table 1, and the numbering of the atoms is shown: 
in Fig. 1. Putting ε--θ and adjusting the resonance integrals, we reduce the dipole moment: 


Table 1. 
a b ο 
A» 1.424À 1.441À 1.445À 
Ryy = Rory 1.377; 1.365А 1.362À 
Rass Ba 1.4364 | 1.4504 | 1.454А 
Rus 1.381 А | 1.365А | 1.361Ä 
σι--αι, 1.073 | 1.047 | 1.032 
gege? | 1.092 | 1.054 | 1.044 
43 1.047 | 1.035 | 0.991 
7 0.622 0.762 0.857 
u 4.8D | 3.0D | 1.9D 


a: Conventional L.C.A.O theory. 
b: L.C.A.O. theory with selfconsistent resonance integrals (e — 0). 


с; L.C.A.O. theory with selfconsistent resonance and Coulomb integrals, &— 1. 


4 


1 


Fig. 1. Numbering of atoms in the fulvene molecule. 
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from 4.8D to 3.0D. This is only about 50% of the correction needed. As the next choice, 
€ was put equal to unity and the resulting L.C.A.O. equations were solved once more by the 
iterational method, starting with the best values of В obtained with e=0 . The convergence 
was poor and the final result was obtained by an interpolation from the first four iterations 
(two overestimates and two underestimates of each charge). The resulting dipole moment, 
1.90), is still too large, which means that ε should exceed 1. 
З In order to estimate the right е, two points were taken into consideration. Firstly, 
from the data compiled by Pullman and Pullman (1952), it appears that, for a wide variety 
of L.C.A.O. orbitals, the total dipole moment of fulvene is a smooth function of gq — the 
charge on the primary atom. On interpolating this function, it is found that the right value 
of the dipole moment (1.2D) corresponds to q4 55 0.928. Secondly, the charges q4 from the 
Once iterated solutions for € = 0, в = 0.7, г = 1.0, and ¢ = 2.0, when plotted as a function 
of e, were found to lie very nearly on a straight line. Assuming that this linear dependence 
holds for further iterations also, we find from the results for в = 0 and ғ = 1 that the correct q4 
and, consequently, the correct dipole moment would be obtained for 


ЕН (4) 


Incidentally, it follows from the plot of q, versus ε that the iterational method diverges for 
values of ε greater than about 1.5. 
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